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Résumé 


La modélisation de la liaison entre le système de fondation d’une structure et son 
site local d’implantation, influence considérablement sur la réponse sismique de la 
structure étudiée. De ce fait, 1l est particulièrement important de considérer l’interaction 
sol-structure dans les zones sismiques où la déformation du sol peut modifier 
radicalement la réponse dynamique de la structure d'un bâtiment. 


Ce phénomène d’interaction sol-structure peut décrire la transmission des mouvements 
sismiques et les effets qui se produisent au niveau de l’interface du contact sol/structure, 
de sorte que, dans l’analyse de la réponse dynamique, on doit inclure les propriétés de la 
structure et également les caractéristiques du sol, et donc, l’interaction sol structure peut 
contribuer à amplifier ou à diminuer la réponse dynamique des éléments de la structure. 


Ainsi, La présente étude a été effectuée en vue d'établir une comparaison entre le 
comportement dynamique des structures en béton armé dont la base est considérée 
comme déformable et celles dont la base est fixe au sol. En premier lieu, on a étudié des 
portiques simples 2D, (4 cas) puis on a fait une étude d’une structure tridimensionnelle. 
Cette comparaison a eu lieu entre les réponses en termes de périodes propres de la 
structure, de déplacements et de l’effort tranchant développé à la base des structures afin 
de mettre en évidence l’importance de la prise en compte de ce phénomène au stade de la 
conception du projet pour différents types de sol 
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Abstract 


The modeling of the bond between the foundation system of a structure and its 
local site of implantation influences considerably the seismic response of the structure 
studied. Thus, it is important to consider the soil-structure interaction in the seismic areas 
where the dynamic response of the soil can significantly change the response of buildings 
structures. 


This phenomena of the soil-structure interaction can describe the transmission of the 
seismic motion and the effects that occur in the plane of the contact soil-structure, so that 
in the analysis of dynamic response , we should include the structure properties and also 
the underlying soil. Therefore, The soil-structure interaction can serve to amplify or 
reduce the dynamic response of the structure. 


So, the present study 1s consisted to compare the dynamic behavior of the Reinforced 
concrete structures with flexible base to that of a fixed base .In the first place, we will 
take an interest in planar frame (4 cases), and then we apply this phenomenon to a 
tridimensional structure. This comparison took place between the responses in terms of 
proper periods of the structure, displacements and the shear stress developed at the base 
of the structures In order to highlight the importance of taking this phenomenon into 
account at the design stage of the project for different types of soil. 
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Introduction générale 


La résolution de tout problème de mécanique nécessite une bonne connaissance de ses 
conditions aux limites. Il est facile d’imaginer, qu’une structure élancée de type poteau 
encastrée à une extrémité, se comporte de manière totalement différente qu’un poteau 
dont l’encastrement est remplacé par une rotule avec ressort. Ses fréquences propres s’en 
trouvent réduites alors que ses déplacements peuvent être amplifiés. Le même 
raisonnement peut être effectué sur un immeuble reposant sur une fondation. Donc, la 
connaissance de la fondation et du sol sous-jacent est essentielle pour bien déterminer le 
comportement réel de la structure. 


Parallèlement, nous pouvons montrer que la présence de la structure influence elle aussi 
le comportement du sol, c’est-à-dire le comportement qu'il aurait en champ libre 
(absence de la structure). En effet, prenons par exemple deux structures, l’une encastrée 
dans un solide rocher, l’autre reposant sur un sol plus souple. La structure encastrée dans 
le rocher ne va pas influencer le champ libre. Le sol étant très rigide, le fait d’y ajouter 
une structure ne change rien à son comportement. La masse de la structure ne change rien 
non plus puisque les efforts d’inerties engendrés par la structure sur la roche rigide n’ont 
pas d’effet sur celui-ci. Par contre, pour la structure reposant sur le sol souple, trois 
phénomènes vont intervenir : le fait d’excaver et d’inclure une fondation rigide va 
modifier le champ libre. Lé fait d'inclure la masse de la structure va induire des efforts 
d’inerties supplémentaires à la base de la structure et donc dans le sol. Ainsi, une 
influence existe bien du sol sur la structure que de la structure sur le sol. C’est les 
problèmes d’interactions entre le sol et la structure. 


En revanche ce phénomène noté généralement ISS, est souvent négligé par les ingénieurs 
malgré son influence sur l’analyse du comportement et le dimensionnement des 
bâtiments. En effet, la répartition des efforts dans les membrures peut être différente suite 
à une analyse qui tient compte de l’interaction sol-structure ce qui affecte sur la sécurité 
du bâtiment. De plus, l’utilisation normale ‘sans ISS” peut être affectée par des problèmes 
de fissurations causées par les tassements différentiels. 


L'objectif de ce travail est d’étudier l’influence dans la prise en compte du phénomène 
d'interaction sol-structure sur la réponse sismique d’un bâtiment en béton armé 


Le travail de ce projet de fin d’étude est présenté en neuf chapitres : 


Y” Le premier chapitre de cette étude est un rappel des phénomènes sismiques et de 
la sismicité au Maroc, 1l comprend quelques séismes historiques au Maroc et des 
généralités sur les séismes pour mieux comprendre le phénomène:;. 

v Dans le deuxième chapitre, on présente le comportement des sols sous un 
chargement dynamique ainsi que l’influence de la nature des sols sur le 
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mouvement sismique et la loi de comportement des sols sous chargement 
cyclique. 

Y” Le troisième chapitre a pour objet de déterminer les caractéristiques dynamiques 
du sol par les différents essais in-situ et dans le laboratoire. 

Y” L’interaction sol-structure fait l’objet du quatrième chapitre; dans ce chapitre on 
montre l’importance de l’interaction sol-structure, la formulation, la modélisation 
d’un problème d’interaction et les méthodes de prise en compte de ce problème. 

Y” Le cinquième chapitre présente les méthodes les plus utilisées pour la le calcul 
des raideurs et d’amortissements des ressorts qui modélisent le sol. 

Y” Le sixième et le septième présentent la modélisation et l’analyse numérique des 
structures en portiques en béton armé, l’étude comporte un calcul dynamique en 
utilisant le code de calcul ETABS qui nous permet de déterminer les périodes, les 
déplacements en tête des portiques et l’effort tranchant à la base des structures. 

vV Le huitième chapitre sera consacré à l’interaction sol-pieu qui consiste à 
modéliser les fondations profondes par des appuis élastiques dont le calcul des 
raideurs sera effectué à l’aide de la du fasciule 62 titre V et l’Eurocode 8 afin de 
comparer les résultats obtenus. 

V Le neuvième chapitre, on va traiter un cas particulier de la liquéfaction dont on va 
parler d’une manière générale de la méthode vérification du risque de la 
liquéfaction. 

v Enfin le travail se termine par une étude de cas sur une structure 
tridimensionnelle. 
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Chapitre I : Sismicité au Maroc 


L Introduction: 


La situation géographique du Maroc, aux frontiéres de la plaque africaine, fait que 
plusieurs régions de notre pays peuvent étre qualifiées comme des zones sismiquement 
actives. Le Maroc est soumis à une activité sismique appréciable à cause de sa situation 
dans un domaine de collision continentale, due à l’interaction entre les plaques 
tectoniques africaine et eurasienne. 


Le Maroc a connu durant son histoire des tremblements de terre destructeurs dont on 
cite en particulier : 


Y” L'effet du séisme de Lisbonne en 1775, tristement célèbre par les désastres qu'il a 
causés, notamment le raz de marrée qu'il a induit et qui a provoqué des 
inondations sur une zone large de 3 Km le long de nos côtes atlantiques ; 

vV Le séisme d’Agadir du 29 Février 1960 de magnitude 5.75 à 6 sur l’échelle de 
Richter et d’intensité X sur l’échelle de Mercalli modifié et qui a provoqué la 
destruction de 70% des constructions et la mort de plus de 12 000 personnes ; 

Y” Le séisme d'Al Hoceima du 26 Mai 1994 de magnitude 5.4 à 5.6 qui a causé la 
mort de deux personnes, de quelques têtes de bétail et la destruction des 
habitations surtout dans le milieu rural (Idsouliyine, Tafensa, Tanda Ifrane, Jbel 
Boubrhel, Arbal et Roudan1) ; 

v le violent séisme d”Al Hoceima du 25 Février 2004 de magnitude 6,5 sur l'échelle 
de Richter, et qui a fait état d'au moins 564 morts et plus de 400 blessés.. 


IL. Généralités sur les séismes 


1. C’est quoi un séisme : 


Le tremblement de terre correspond à une vibration du sol provoqué par une 
libération soudaine de l’énergie de déformation accumulée dans la croûte terrestre ou 
dans la couche sous-jacente appelé manteau. 


Lorsque les contraintes dépassent un certain seuil, une rupture d’équilibre se produit et 
donne naissance aux ondes sismiques qui se propagent dans toutes les directions et 
atteignent la surface du sol. Ce mouvement du sol excitent les ouvrages par déplacement 
de leurs appuis et sont plus ou moins amplifiés dans la structure. Le niveau 
d'amplification dépend essentiellement de la période de la structure et de la nature du sol. 


2. Terminologie sismique : 
2.1.Foyer et épicentre : 
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Les ondes sismiques semblent provenir d’une région déterminée appelée foyer ou 
hypocentre du séisme. Le point de la surface du sol situé à la verticale du foyer est son 
épicentre. La distance entre ce point et l’endroit d’enregistrement de l’accélérogramme 
est appelée distance épicentrale. 
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Figure 1: position du foyer, l'épicentre et la faille 


La détermination de ces positions est un probléme complexe et donne des résultats plus 
au moins précis car la libération de l’énergie n’est pas un évènement ni ponctuel, ni 
instantané. 


2.2.Sismographe et sismogramme : 


Le sismographe est l’appareil qui mesure le mouvement du sol et l'enregistre sur un 
support visuel. Le tracé de ce mouvement s'appelle un sismogramme. 





Figure 2: sismographe 


3. Les caractéristiques des séismes : 
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Le séisme est caractérisé par : 
3.1.Son intensité : 


Elle mesure les effets et dommages du séisme en un lieu donné. Ce n’est pas une mesure 
objective, mais une appréciation de la manière dont le séisme se traduit en surface. On 
utilise habituellement l’échelle MSK, qui comporte douze degrés. 





, Effets du tremblement de terre sur 
Degré Secousse "FP 
personnes 
imperceptible | nonressent | 


ressenti 
rarement : 
. sporadiquement 
perceptible 2 
(étages supérieurs) 





ressenti par quelques À 
| par qued léger balancement des objets 
faible personnes au repos NT A lusiees 
dans leur habitation P 
som qe o 


maisons ; quelques 
dormeurs éveillés 


largement 
observée 


vibration des portes, fenêtres et 
assiettes 


beaucoup de bâtiments crevasses dans 
IX destructions panique s'effondrent selon la qualité de la | terrain : beaucoup 
construction de glissements 


beaucoup de constructions rallais tordus, 
très destructive | panique générale s'effondrent quelle que soit la | glissement de terrain 
qualité sur versants 


, l la plupart des constructions modifications du sol 
XI dévastatrice , À 
s'effondrent : tsunami 


pratiquement toutes les structures 


bouleversement de 
au-dessus et au dessous du sol 


XII | catastrophique 


la topographie : 


sont fortement endommagées ou . 
tsunami 


détruites 
Tableau 1: les intensités selon l'échelle de MSK 
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Différentes échelles d’intensité ont été définies, quí classent les effets sismiques suivant 
leur importance croissante, à l’aide de descriptions conventionnelles : 

v L’échelle macrosismique internationale d’intensité, et l’échelle Mercalli. 

v L’échelle macrosismique MSK. 

v L’échelle MM (Mercalli Modifiée) utilisée aux États-Unis. 

v L’échelle Japonaise JMA (Japan Météorologique Agency). 

v L’échelle macrosismique Européenne EMS. 


3.2.8a magnitude : 


Elle traduit l’énergie libérée par le séisme. Elle est généralement mesurée par 
l’échelle ouverte de Richter. 


Magnitude Effets engendrés 


Destruction totale à l’épicentre, et possible 
sur plusieurs milliers de km 
Dégâts majeurs à l’épicentre et sur 
plusieurs centaines de km 
7 Importants dégâts à l’épicentre, secousse 
ressentie à plusieurs centaines de km 
Dégâts à l’épicentre dont l’ampleur dépend 
de la qualité des constructions 


dommages mineurs près de l’épicentre 


Secousse sensible, mais pas de dégâts 


sensible pour la plupart des gens 
Secousse ressentie uniquement par des 
gens au repos 


Secousse imperceptible 


Tableau 2: les magnitudes selon l'échelle de Richter 





Remarque : Relation entre intensité et magnitude 


Il n’y a pas de relation directe entre l’intensité et la magnitude. Les deux grandeurs sont 
difficilement comparables. Un séisme de forte magnitude avec un foyer profond et dans 
une région peu peuplée sera peu destructeur et donc sera qualifié de faible intensité. Au 
contraire, un séisme superficiel, même de magnitude moindre pourra être très destructeur 
et donc caractérisé par une grande intensité 


3.3.La fréquence et la durée des vibrations : 


Ces deux paramètres ont une incidence fondamentale sur les effets en surface. 
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3.4.L”accélérogramme: 
L’accélérogramme représente l’évolution en fonction du temps de l’accélération du sol. 
3.5.Le spectre de réponse: 


La pratique actuelle la plus répandue consiste à définir le chargement sismique par un 
spectre de réponse.les spectres caractérisent les séismes beaucoup plus clairement que les 
accélérogrammes. 

L'utilisation d’un spectre de réponse s'effectue différemment selon qu'il s’agit d’un 
spectre de déplacement ou d’un spectre d’accélération. La définition de l’action sismique 
par un spectre de réponse est une méthode particulièrement adapté pour 


> synthétiser par une enveloppe couvrante tout le domaine des fréquences les efforts 
de plusieurs accélérogammes aux contenus très différents ; 

> Prendre en compte la nature du terrain situé sous la construction ; 

> Tenir compte de l’intensité probable du séisme, du niveau de protection recherché 
et ceci en fonction de l’importance du bâtiment étudié, ce quí revient à caler le 
spectre au niveau de l’action sismique. 


3.6.Aléa sismique : 


L’aléa est la manifestation d’un phénomène naturel d’occurrence et d’intensité 
donnée. Appliquée à la sismicité, l’analyse de l’aléa étudie l’occurrence des 
tremblements de terre et les accélérations ou mouvements du sol qui en découlent. 


3.7.Le risque sismique : 


C’est la conséquence de l’aléa sismique dont les effets peuvent mettre en jeu un 
grand nombre de personnes et occasionner des dégâts importants. 


4. Propagation des ondes sismiques : 


Un tremblement de terre correspond à une vibration transitoire du sol provoquée par les 
ondes sismiques émises lors d’un déplacement de sol (rupture) le long d’une faille active, 
à partir du foyer. 


Les vibrations lors d'un séisme se propagent dans toutes les directions. On distingue deux 
types d'ondes, les ondes de volume qui traversent la Terre et les ondes de surface qui se 
propagent à sa surface. Sur les enregistrements des sismographes, elles se succèdent ou se 
superposent. Leur vitesse de propagation et leur amplitude sont modifiées par les 
structures géologiques qu'elles traversent, c'est pourquoi, les signaux enregistrés sont la 
combinaison d'effets liés à la source, aux milieux traversés et aux instruments de mesure. 
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4.1.Ondes de volume: 


Elles prennent naissance dans le foyer et se propagent à l’intérieure de la terre sous 
deux formes: 


4.1.1. Ondes primaires ou ondes longitudinales (P): 


Ce type d'ondes fait osciller les ouvrages d'avant en arriére dans la méme direction que 
leur propagation, elles induisent donc, un mouvement vertical dans les ouvrages. Ce 
mouvement sera d'autant plus important que l'ouvrage se trouve a proximité de 
l’épicentre ; Elles se propagent dans la croûte terrestre avec une vitesse de 6 à 8km/s et 
s’accompagnent d’un changement de volume (compression et dilatation alternées). 


Dilatations Compressions Milieu normal 


a == == 


PER E 





Figure 3: Onde de compression (P) 


4.1.2. Onde secondaire ou ondes de cisaillements (S): 


Elles se propagent dans la croûte terrestre, avec une vitesse d'environ 3 a 5km/s et 
s'accompagnent d’une distorsion dans le plan perpendiculaire à la direction de 
propagation, provoquant un cisaillement sans changement de volume. C’est bien 
ce type d’ondes qui produit un mouvement horizontal dans les ouvrages. 


On peut calculer la vitesse de propagation de ces ondes à partir des 
caractéristiques du sol: 


Ondes 5 





Figure 4: Onde de cisaillement (S) 
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4.2.Ondes de surface 


Les ondes de volume qui arrivent à la surface de la terre donnant naissance à des ondes 
de surface de deux types: 


4.2.1. Ondes de Rayleigh(R) 


Ce sont des ondes pour lesquelles les points du sol décrivent des ellipses dans le plan 
vertical et dans le sens rétrograde. Ce mouvement entraîne des compressions ou des 
tractions, ainsi que des cisaillements dans le sol. 


AAA 


e 

PRO 
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Ondes de Rayleigh ————— a 





Figure 5: Onde de Rayleigh (R) 
4.2.2. Ondes de Love (L) 


Ce sont des ondes pour lesquelles les points du sol se déplacent dans un plan 
tangent à la surface et perpendiculairement à la propagation. Leur mouvement 
n’engendre que des contraintes de cisaillement. 





Ondes de Love sm 


Figure 6: Ondes de Love (L) 


Les ondes de cisaillement (S) et (L) sont plus dangereuses pour les constructions 
que les ondes (P). 
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II. Zonage sismique: 


v [Pour simplifier le calcul des charges sismiques et uniformiser les exigences de 
dimensionnement des structures à travers de grandes régions du pays, le RPS 2011 
utilise l'approche des zones. Il s’agit de diviser le pays en plusieurs zones de sismicité 
homogène et présentant approximativement le même niveau d’aléa sismique pour une 
probabilité d’apparition donnée. 

Y” Dans chaque zone, les paramètres définissant l’aléa sismique, tels que l’accélération 
ou la vitesse maximale horizontale du sol, sont considérées constantes. 

Y” Pour identifier adéquatement le caractère particulier d’un séisme en un endroit donné 
le règlement adopte un zonage séparé pour les paramètres Amax et Vmax .exprimées 
respectivement en fraction de 1g et 1 m/s. 

Y” Chacune des deux cartes de zonage sismiques adoptées par le RPS 2011 comporte 
actuellement cinq zones (0 à 4) reliées à l’accélération horizontale maximale du sol 
Amax (Za=0;1;2;3;4)etaà la vitesse horizontale maximale du sol Vmax (Zv = 0 ; 
1,2 ;3 ;4), pour une probabilité d'apparition de 10% en 50 ans , ce qui correspond à 
une période de retour de 475. Cette probabilité est considérée raisonnable, car elle 
correspond à des séismes modérés, susceptibles de se produire plusieurs fois dans la 
durée de vie d’une structure. Les deux cartes des zonages sismiques relatives à Amax et 
à Vmax Sont présentes respectivement dans les figure 7 et 8 respectivement. 


Vitesse (cm/s) 





Figure 7: Zonage sismique en vitesse pour des probabilités de 10% en 50 ans Maroc 
2011 (vitesse cm/s) 
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Figure 8: Zonage sismique en accélération pour des probabilités de 10% en 50 ans 
Maroc 2011 (accélération %g) 
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Chapitre II : Phénomène de liquéfaction 


L Explication du phénomène : 


Les sables fins saturés en eau peuvent être le siège d'un phénomène encore plus 

spécifique au contexte sismique : la liquéfaction des sols. L*explication du phénomène 

est la suivante : 

La résistance au cisaillement d'un sol est donnée par la relation : 
T=C+(0—U)Xtan(qp) 

Avec: 

c: Cohésion du sol. 

u : Pression interstitielle dans le sol. 

p: L'angle du frottement interne du sol. 


Lorsqu'un séisme frappe une région, 1l donne lieu a des vibrations qui vibrent le sol, ces 
vibrations réarrangent les grains du sol entre eux. L’eau emprisonnée entre les grains veut 
sortir suite à ce réarrangement, mais comme la durée du séisme est faible, cette eau ne 
trouve pas le temps nécessaire pour s’échapper, ce qui conduira à l’augmentation de la 
pression interstitielle. 


Il s’agit d’un phénomène géologique qui se produit pour un sol pulvérulent saturé (sables 
et limons et argiles sensibles) quand une perte soudaine de force et de rigidité est vécue, 
entraînant parfois des grands déplacements permanents du sol et permettant ainsi 
l'enfoncement des objets lourds situés en surface : 


Ce phénomène se produit pendant le mouvement sismique et engendre : 


e la diminution de la rigidité du sol, 

e et/ou la diminution de sa résistance au cisaillement, due à 
l’augmentation de la pression de l’eau interstitielle, 

e des déformations permanentes significatives (tassements, 
glissements), 

e voire une quasi-annulation des contraintes effectives (étalements). 


Comme exemple de sol liquéfiable on trouve les sables fins et les limons, peu compacts 
et saturés, ou les surpressions interstitielles engendrées par les sollicitations sismiques 
peuvent conduire à la perte de consistance du sol, phénomène appelé liquéfaction. Des 
glissements de terrain, des enfoncements et des renversements d’immeubles en sont les 
conséquences visibles 
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II. Justification à la liquéfaction 
L’approche utilisée pour l’évaluation de potentiel de liquéfaction d’un sol est celle 
proposée par Seed et Idriss on suivant les étapes suivantes : 


1. Identification de couches potentiellement liquéfiables 


Les paramètres déterminant la liquéfaction des sols sont : 
e La granulométrie. 
e La forme des grains. 
e Le poids volumique du sol en place. 
e La contrainte effective, due essentiellement au poids propre du sol. 


Selon le RPS 2011, les sols susceptibles à la liquéfaction sont : 


= Les sables et les limons : 
- Un degré de saturation S, ~ 100% 
- Une granulométrie caractérisée par: 


Un coefficient d'uniformité Cy < 15 . 
Un diamètre 0.05mm < Ds < 1.5mm . 


= Les sols argileux fins : 
- Un diamètre D:, > 0.005m 
- Une limite de liquidité LL < 35% . 
- Une teneur en eau naturelle Wp > O.9LL . 
- Un indice de liquidité <0.75. 
= Les sols sableux dont la courbe granulométrique s’inscrit dans le fuseau des sols 
liquéfiables. 


La classification des sols peut se faire par des méthodes géophysiques ou des méthodes 
géotechniques : 


1.1.Méthodes géophysiques : 


La classification dans ce cas nécessite la connaissance de la vitesse des ondes de 
cisaillement Vs. La vitesse Vs est calculée à l’aide d’un essai qui consiste à assimiler 
l’onde séismique traversant un sol donné par un coup de marteau, en fait, on mesure le 
temps durant lequel l’onde traverse une distance x et on calcule la vitesse par 
l’expression : 
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1.2.Méthodes géotechniques : 
La classification dans ce cas nécessite la connaissance de l’un des paramètres suivant : 


- La résistance à la pénétration standard N60 calculée à l’aide de l’essai SPT : 
- La résistance du sol non drainé au cisaillement Su déterminée par l’essai triaxial 


Lorsque l’un de ces paramètres est calculé, la classification se fait alors à l’aide du 
tableau suivant : 


CLASSIFICATION DES SITES 

Résistance du Vitesse des 
sol non drainé ondes de 

au cisaillement | cisaillement Ys 


Résistance à la 
TYPE DE SOL Pénétration 
Standard No 


Rocher sain à moins de 3 m des 
fondations 


Rocher altéré; sols cohérents 
trés raides, sols pulvérulents 
très denses, marnes ou argiles 
très consolidés 

Sables et graviers 
moyennement compacts, argile 
moyennement raides 

Sols pulvérulents láches 


Tout Sol de profondeur de plus 
de 3 m ayant Teneur en eau 
W> 40% ; Indice de plasticité 
PI > 20 


Sols à conditions 
exceptionnelles nécessitant une 
étude appropriée de leurs 
caractéristiques : 

1. Les argiles ayant un PT > 75 
et plus de $ m d'épaisseur ; 

2. Les argiles raides de molles à 
moyennes dont l'épaisseur est 
supérieure 430 m ; Mécessite des études géotechniques 
3. Les sols susceptibles d'être 
liquéfiables ; 

4. Les sols susceptibles 
d'affaissement sous des 
sollicitations sismiques ; 

5. Les sols avec grande teneur 
en matière organique sur une 
épaisseur de plus de 3 métres. 





Tableau 3: Classification des sites (selon RPS2011) 
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2. Calcul du coefficient de sécurité contre la liquéfaction du site : 


Le coefficient de sécurité est calculé par la formule suivante : 





E CRR 
S = (CSR 
Avec : 
CSR : Ratio de résistance cyclique 
T 
CSR = — 


O'y0 
Tay : Contrainte de cisaillement moyenne. 
O'yo : Contrainte effective à une profondeur donnée. 
CRR : Rapport de résistance cyclique. 


2.1.Evaluation de la contrainte de cisaillement moyenne : 


Il est possible d’utiliser la méthode simplifié de Seed - Driss, 1971 pour la détermination de la 
contrainte de cisaillement maximale : 


Avec : 

(Tmax)r : Contrainte de cisaillement maximale due au séisme à cette profondeur ; 
Oyo : Contrainte verticale totale à la profondeur considérée; 

Amax : Accélération maximale du sol en surface; 

g : Accélération de la pesanteur. 


Or en réalité la contrainte de cisaillement maximale sera inférieure à celle calculée par 
l’équation c1-dessus si on prend en considération la flexibilité du sol. 


(Ja = Ta X (Taueli 
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Figure 9: procédure de la détermination de la contrainte maximale de 
cisaillement 


Où r4 est un coefficient réducteur de contrainte de cisaillement, ce coefficient égale 1 en 
surface et décroit avec la profondeur. 


Seed & Driss ont fourni une formule approximative qui exprime la variation de rg : 
Ya = 1 — 0.00765zZ Pour z<915m 
Ya = 1.174 — 0.0267z Pour 9.15m < z < 23m 


Or l'équation trouvée est pour les sollicitations sismiques d’amplitude constante, donc on 
doit nécessairement convertir le diagramme T(t) en un diagramme « équivalent » pour N 
cycles. Cependant après avoir effectué plusieurs exemples, ils ont trouvé avec un degré 
de sécurité raisonnable la contrainte de cisaillement moyenne équivalente est environ de 
65% de la contrainte de cisaillement maximale 


Tav = 0.65 X (Tmax)a 
2.2.Evaluation du rapport de résistance cyclique : 


La meilleure approche pour estimer le rapport de résistance sismique CRR est d'extraire 
des spécimens de sol afin de les faire tester dans le laboratoire. Cependant, une telle 
approche se heurte au problème de représentativité de l'échantillon par rapport à l'état 
initial du matériau en place, soit pour son état de contrainte initiale qui est difficilement 
atteint au laboratoire, soit parce que l'échantillon est assez remanié par les techniques de 
forage et d'échantillonnage, pour que les résultats de l'expérimentation soient 
significatifs. 


En utilisant, exclusivement, quelques techniques spécialisées de prélèvement telles que la 
congélation du sol, on peut tirer des spécimens assez représentables pour les essais au 
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laboratoire. Hors, le coût abusif de telle technique a limité son application dans les projets 
de grande importance. 


Pour éviter tout problème associé au prélèvement et aux essais de laboratoire, des essais 
en place sont devenus, de plus en plus, l'état de pratique d'une investigation courante de la 
liquéfaction. Outre le SPT (Standard Pénétration Test), on peut citer le CPT (Cône 
Pénétration Test) et l'essai géophysique en onde Vs (vitesse d'onde de cisaillement).et le 
tableau suivant explose les avantages et les inconvénients de ces différents essais : 


r a Type d'essais 
Particularité 


Mesures effectuées auparavant 
i e | nombreux nombreux nombreux 
dans les sites de liquéfaction 


Partiellement 
drainé, grande 
déformation 


Type de comportement contrainte-déformation 
qui influence sur l'essai 


Drainé, grande Petite 
déformation déformation 


e a oa Médiocre a | 


Bonne pour les 

Détection de la variabilité des essais , 

a nooo. Très bonne Assez bonne 

dépôts de sol. étroitement 
espacés 
Types de sol dans lesquels l'essai 
YP a | Non graveleux | Non. Graveleux | Tous les types 
et recommandés 


Échantillon de sol récupéré | oui | Non | non | 


Tableau 4 : Avantages et inconvénients des différents essais d'estimation en place E 
la résistance a la liquéfaction 





3. Evaluation de CRR par Pessai de Pénétration au carottier (SPT) : 


* Méthode graphique de Seed: 
CRR = CRR> 5 X MSF 
Avec : 
MSF : Facteurs d'échelle d’ampleur MSF 


Tenant compte de la magnitude du séisme et donné par le diagramme 
suivant : 
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Figure 10: Les corrections apportées au coeffcient MSF par des différents 
investigateurs 


Les essais d'évaluation de la résistance à la liquéfaction basés sur les résultats du SPT ont 
été assez robustes pendant des années. Ces essais sont introduits en grande partie dans le 
diagramme de Seed par les couples (N1)60- CRR représenté dans la figure suivante : 


FINES CONTENT 2 5% 


05846 > Modified Chinese Code Proposal (clay content = 5%)@ 


Cyclic Stress Ratio (CSR) or Cyclic Resistance Ratio (CRR) 


Liquefaction penis Liquefaction 


0” A 
g | ¡Adjustment Pan - America data m 
Recommended]| Japanese data e o 
By Workshop Chinese data A 


30 
Corrected Blow Count, (Ny)¿g 





Figure 11: Le rapport de résistance cyclique en fonction de nombre de coups corrigé et 
de la teneur en fine 
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A L'université du Texas A. F. Rauch (communication personnelle, 1998) a approché la 
courbe de base du sable propre indiqué dans la figure par Ass suivante : 


z 1 (N1)60 50 
CRR;55 = na a 00 Valable pour (N; )¿y < 30 


Ils existent plusieurs paramètres influant la valeur du nombre N de coups SPT : 

e La structure du matériau 

e Son degré de surconsolidation 

e Sa densité relative et la contrainte verticale effective qu'il 

supporte 

Une correction approximative de l'influence des fines sur le CRR doit être faite. Ces 
équations ont été développées par I.M. Idriss avec l'assistance de R.B. Seed pour corriger 
(N1)60 à une valeur équivalente (N1)60 du sable propre en tenant compte de la teneur en 
fine : 


(Ni)socs = æ + P X (N:)60 


a et p sont déterminés à partir des relations suivantes : 
a=0Etf=1 pour FC < 5% 


a = exp(1.76 — 7) et B = 0.99 + FC! pour 5% < FC < 35% 


a =5 Etf =1.2 pour FC > 35% 
Outre que la teneur en fines, plusieurs facteurs influencent les résultats du SPT, le tableau 


suivant ainsi que l'équation introduit les corrections suivantes: 
(Ni)so = Nm X Cy X Cp X Cp X Cp X Cs 


Nm : résistance à la pénétration mesurée. 

CN: facteur de correction qui normalise Nm à une contrainte de confinement de100 KPa. 
CE : correction pour l'énergie du marteau. 

CB: facteur correcteur du diamètre du trou de sondage. 

CR : facteur de correction sur la longueur de la tige du forage. 

Cs : correction sur la procédure de prélèvement. 

Ou facteur correcteur sur la tarière (avec ou sans cale) 


Variables de l'équipement 


Paragraphe 
Pression de confinement srap 
ao 


Presion de confinement | | q | <17 


n : Marteau avec tore 
Diamé 
sondage 
Diamé 
sondage 
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Diamètre du trou de 


200 mm Cs 1.15 
sondage 


<3m | 
longueur de la tige | 34m | CR | 08 | 


Longueur de la tige 10-30m 


Méthode 
, Echantill 1. 
d'échantillonnage d 
Méthode 
| | Echantill li 1.1-1. 
d'échantillonnage AAA 


Tableau 5: Corrections sur le SPT (modifiées par skemptin 1986) comme présentées par 
Robertson et Wride (1998) 





Comme la résistance à la pénétration (SPT) augmente en augmentant la pression de 
confinement, un facteur de correction de pression est appliqué (Seed et Idriss 1982). Ce 
facteur est généralement calculé par l'équation: 
la 0.8 
Cy = C) 
O vo 
CN ne doit pas excéder la valeur de 1.7, Kayen et al (1992) ont suggéré l'équation 
suivante, qui limite la valeur maximale de Cy à 1.7, et selon leur opinion elle fournit une 
meilleure conformité à la courbe spécifiée par Seed et Idriss (1982) : 
g 
Cy = 2.2/(1.2 +35) 


a 








Avec : 
O'y0: Contrainte effective pour une 
P,: Pression atmosphérique (100Kpa) 
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Figure 12: Facteur correcteur CN de la contrainte effective en fonction du nombre de 
coups N-SPT 
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Chapitre III : Comportement des sols sous chargement 
cyclique 


L Introduction 

La loi de comportement est un moyen indispensable dans l’étude des problèmes 
d'interaction sol structure, c’est un modèle théorique décrivant le comportement du 
sol et permettant de relier le tenseur de déformation € au tenseur de contrainte ©. 
La modélisation du sol nécessite le choix judicieux d’une loi de comportement 
réaliste, approchant mieux le comportement du milieu .Le choix du modèle approprié 
dans l’analyse de l’interaction sol-structure n’est pas unique et dépends de plusieurs 
facteurs : 


Y” Type de sol 
Etat de sol 
type de la structure en contact avec le sol 


SN SS 


nature des charges transmises 


Le modèle ainsi déterminé permet d’évaluer le comportement du sol soumis à des 
chemins de contraintes quelconques, tels que ceux suivis in-situ lors d’un séisme. Le 
modèle doit être formulé en contraintes effectives puisque le comportement du sol est 
régi par ces contraintes. 


élastique-linéaire elastique-non linéaire rigide-par fa ilemenl 
plastique 





élasto-plastique élasto-plastique élasto-plastique fragile Hyper bolique 
sans écrouissage avec ÉCTOUISSALE (Kondner-Duncan) 


© Critère de rupture - loi de plastification 
(par ex. : Mobr-Coulomb, Drucker-Prager, von Mises) 





Figure 13: Lois de comportement 
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II. Chargement monotone 


On doit tout d’abord rappeler le comportement des sols sous chargement quasi-statique 
monotone pour mieux comprendre son comportement chargement dynamique. 

Les courbes efforts déformation obtenues à l’appareil triaxial pour trois chemins de 
contraintes différents: 


% Compression isotrope; 

% Essai à déformation latérale nulle: 

% Essai triaxial classique à contrainte radiale (ou pression de confinement) 
constante. 


L”examen des résultats expérimentaux montre que: 


Il existe un domaine de faible déformation pour lequel la relation entre la 

contrainte appliquée et la déformation associée est linéaire. 

e Au-delà d’un certain seuil de contrainte et quel que soit le chemin de contrainte 
suivi, la relation effort- déformation n’est plus une droite. Le comportement du 
sol cesse d’être linéaire. Il faut noter que le non linéarité n’exclurait pas à priori 
Phypothèse d’un comportement élastique qui se traduirait par un trajet de 
décharge identique à celui de la charge. 

e Si à partir d'une contrainte suffisamment élevée on effectue une décharge de 
échantillon, le trajet suivi lors de la décharge n’est plus identique au trajet suivi 
lors de la charge. Il y apparition d’une déformation résiduelle. 

e Pour certains trajets de chargement la déformation devient très importante pour 

une valeur finie de la contrainte appliquée. Le sol atteint un état de rupture. 


Sous chargement quasi-statique monotone le comportement du sol est donc caractérisé 
par un domaine d'élasticité (domaine à l’intérieur duquel les déformations restent 
élastiques) 


IHI. Chargement cyclique 


1. Le comportement du sol sous chargement cyclique : 


Le substratum, plan horizontal, se comporte comme une table vibrante dont le 
mouvement en surface doit coincider avec le mouvement sismique donné. Etant situé en 
profondeur, son mouvement est très différent de celui en surface du sol (en champ libre) 
et n’est pas connu a priori. 


En réalité, Les contraintes obtenues de l’essai triaxial de compression ne simulent pas 
parfaitement la condition du chargement : 
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Figure 14: Transmission des ondes du cisaillement dans une couche du sol depuis le 
susbtratum 


Il est d'usage dans les calculs de réponse dynamique d'un profil de sol, ou dans les 
problèmes d'interaction sol-structure, de considérer que le mouvement sismique a pour 
origine une onde de cisaillement se propageant verticalement depuis le substratum. Dans 
ces conditions, on prend un élément de sol initialement soumis au poids des terres au- 
dessus : 





Figure 15: Séquence de chargement idéalisée 


Initialement, dans le cas d'un profil de sol horizontal, l'élément est en équilibre sous les 


contraintes verticales effectives On” et horizontales effectives Ko Oo’, où Ko est le 
coefficient de restitution ou coefficient des terres au repos. 

Lorsque ce coefficient n’est pas déterminé expérimentalement, on peut considérer les 
formules suivantes : 


Y 


% Dans un milieu élastique qui aurait été chargé verticalement sans déformation 


latérale : Ko = — 
Tel que v est le coefficient de poisson du milieu. 
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4 


* Dans un terrain argileux, normalement consolidé, on peut appliquer la formule de 
JAKY : Ko = 1 — sin(ọ) 


Le passage de l'onde de cisaillement se traduit par l'application sur les faces horizontales 
de l'élément de sol, et donc sur les faces verticales pour maintenir les conditions 
d'équilibre, d'une contrainte de cisaillement T(t), sous l'effet de cette contrainte, 
l'échantillon subit une déformation de cisaillement simple, qui pour un matériau à 
comportement élastique, se traduirait par une variation de volume nulle. La déformation 
de cisaillement, appelée distorsion, est définie par : 


õu Au 
Y az Ah 








Figure 16: Séquence de chargement idéalisée 


L'enregistrement d’un courbe effort — déformation T = f (y) montre que pour un cycle 
fermé (figure 17) le comportement du sol, est caractérisé par une boucle appelée boucle 
d’hystérésis, dont la surface et l’inclinaison dépendent de l’amplitude de la déformation 
au cours du cycle. 


OL us 


Courbe de 
1°" chargement. 


Contrainte T 





Distorsion Y » 





Figure 17: Courbe effort déformation cyclique 


Il est commode et classique de définir cette boucle à l’aide de deux paramètres: 


v La pente de la droite joignant les deux extrémités de la boucle est le 
module sécant: Gs 
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v Parede la boucle est l’énergie dissipée qui caractérise le coefficient 
d'amortissement du milieu. 

v La partie hachurée désigne l’énergie dissipée si le comportement est 
élastique 


Dés que le chargement cyclique n’est pas fermé, le comportement devient plus complexe 
à décrire. Un exemple en est donné sur la figure 14, Jusqu'en B, le trajet suivi est 
identique à celui de la figure 13 (courbe de premier chargement puis courbe de décharge). 
Au point B, tel que Op < Oa, le signe du chargement est de nouveau inversé; le trajet 
suivi est donné par la courbe BC puis éventuellement par CE si le chargement change à 
nouveau de signe en C. Si au contraire, le chargement est poursuivi au-delà de C, le trajet 
est représenté par CA puis au-delà suit de nouveau la courbe du premier chargement. 





Figure 18: chargement cyclique quelconque 


La courbe représentée dans la figure est centrée, en réalité elle a la forme suivante : 








Figure 19: Chargement cyclique non centré à l’origine 
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2. Les domaines de comportement 


Le sol étant constitué d'un squelette granulaire, les déformations d’un échantillon 
peuvent provenir à la fois d’une déformation des grains et de leur réarrangement. Il peut 
se comporter de façon très différente : tantôt comme un corps élastique, tantôt comme un 
matériau élasto-plastique ou plastique parfait. Plusieurs critères sont proposés pour 
évaluer le comportement des sols définis à partir du type de charge, du domaine de 
déformation. En particulier les déformations de cisaillement susceptible de se développer 
lors d’un séisme. Le tableau 6 illustre les différents domaines de comportement ainsi que 
la description des modèles mathématiques qui peuvent être utilisés dans l'analyse 
numérique. 


non a Inférieur à 10° Entre 10° et 10* Supérieure à 10° 
cisaillement y 


Loi de comportement Hystérétique stabilisé | Fortement non linéaire 


Tableau 6: Modèles de comportement et méthode d'analyse en fonction de 
l'amplitude de la déformation cyclique (d’après Davidovici 1985) 





On classifie le comportement des sols en quatre domaines distincts de comportement 
selon la valeur de l’amplitude de déformation pour n’importe quel type de sollicitation 
cyclique : 


2.1.Domaine linéaire (0 < ye < 107) 


Pour des cycles d’amplitude inférieure à 10° autour d'un état contrainte-déformation 
donné, le comportement du sol peut être considéré comme élastique linéaire : aucun 
réarrangement irréversible ne se produit. Le comportement peut être considéré comme 
élastique quasiment linéaire. Toutefois, le comportement du sol n’est pas totalement 
élastique linéaire en raison d’un amortissement non nul. Cependant la valeur du 
coefficient d’amortissement est très faible. 


2.2.Domaine hystérétique stabilisé (10° < ye < 107?) 


Le domaine des petites déformations ou domaine hystérétique stabilisé (10 °< ye <10*), 
où le comportement du sol est hystérétique jusqu’à un niveau de déformation 10*. Les 
boucles des cycles contrainte-déformation sont en effet, nettement plus ouvertes. 
Cependant ces boucles se stabilisent, c'est-à-dire qu’elles adoptent la même forme quel 
que soit le nombre de cycles. Les modules d’Young et de cisaillement décroissent avec le 
niveau de déformation. L’amortissement augmente avec la déformation et peut 
atteindre 7%. 
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On peut adopter pour ce domaine un comportement élasto-plastique qui présente des 
déformations réversibles (élastiques) et aussi des déformations irréversibles (plastique). 


2.3.Domaine hystérétique non stabilisée (104 < ye < 103): 


Quand l’amplitude de déformation dépasse 10*, le phénomène non-linéaire devient de 
plus en plus important. On constate alors une accumulation de déformation volumique 
(un réarrangement des grains (rotations, glissements) selon le chargement. Ce 
réarrangement risque d’induire une anisotropie dans le matériau) en conditions drainées 
et une augmentation de la pression interstitielle en conditions non drainées. 
L'augmentation de la pression interstitielle peut alors mener à la liquéfaction, un 


phénomène très important en mécanique des sols au cours du quel le sol perd sa 
résistance de cisaillement. 


2.4.Domaine des lois incrémentales (domaines de grandes déformations) 
(Ye > 10°) 


Enfin, dans ce domaine, la sollicitation est susceptible de conduire à la rupture des 
grains, et donc modifier profondément les caractéristiques du matériau. De plus, 
l’irréversibilité inélastique et l’effet visqueux deviennent prédominants. L’effet visqueux 
caractérisé par la vitesse de sollicitation peut jouer un rôle particulièrement pour les 
argiles molles. Seule l’approche incrémentale associée avec l’intégration étape par étape 
le long du chemin de sollicitation permet de décrire ces phénomènes variés, qui relie à 
chaque instant le tenseur complet des taux de contraintes effectives au tenseur des taux de 
déformations. Le coefficient d'amortissement se stabilise et tend vers une valeur 
maximale. Les modules d’Young et de cisaillement sont très faibles par rapport à ceux du 
domaine des très petites déformations. 


La courbe (figure 20) illustre les différents domaines de comportement en fonction 
de l’indice de plasticité et la distorsion. 
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Figure 20: Seuils de déformation délimitant les domaines de comportement 
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3. Le comportement adopté : 


Pour notre étude on va s’intéresser au comportement élastique-linéaire du sol : 





Figure 21: comportement élastique-linéaire 


Ce modèle n’est valable pour le sol que dans un domaine très restreint et pour 
des contraintes très faible. Le comportement élastique des sols est limité seulement 
aux très petites déformations (y < 107). De plus, ce modèle s’accorde avec notre étude en 
raison que notre pays est de faible sismicité et donc, les contraintes sont très faibles. 


En générale une loi élastique a la forme suivante : 


e [|£] est le tenseur de déformation 

e [|o] est le tenseur de contraintes 

e [C]est un tenseur du 4°" ordre contenant 36 composantes distinctes dans le cas 
anisotrope mais dans le cas d’isotrope, 2 composantes sont suffisantes. 


Dans le cas unidirectionnel et pour des contraintes de cisaillement cette loi a la forme : 
T = G.y 


e Tt:Contrainte de cisaillement 
e G: module de cisaillement 
e y: déformation de cisaillement (distorsion) 


Conclusion 

La description réelle du comportement des sols sous chargement cyclique nécessite la 
détermination d’une loi de comportement globale qui prend en compte le caractère non 
linéaire et les déformations irréversibles du sol. Et comme dans notre étude, le sol subi à 
des faibles déformations en raison que notre pays est de faible sismicité d’une vision 
générale, on pourra donc adopter une loi de comportement élastique qui permet de 
caractériser le sol avec un minimum de paramètres. 
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Chapitre IV : Les caractéristiques dynamiques du sol 


L Introduction : 


Le problème des mesures des caractéristiques des sols constitue l’un des aspects 
fondamentaux de la mécanique des sols, et plus spécifiquement, de la dynamique des 
sols. Les modèles les plus élaborés, les calculs les plus compliqués ne sont d’aucune 
utilités si les paramètres à entrer dans ces modèles, ou ces calculs, sont incorrects ou mal 
connus. Par hypothése, le sol est supposé étre un milieu semi-infini, homogéne et 
isotrope, de comportement élastique et linéaire, les couches du sol étant supposées 
horizontales. Le sol est défini par ses caractéristiques mécaniques : 


% le module dynamique d* Young E, ou le module transversal G dynamique ; 

% le coefficient de Poisson dynamique ; 

% la masse volumique ; 

% l'amortissement interne de frottement ; 

Généralement, la connaissance des caractéristiques de sol est assez imprécise. 
L'hypothèse faite sur l’homogénéité est rarement vérifiée : les sols sont la plupart du 
temps stratifiés et leurs propriétés mécaniques augmentent généralement avec la 
profondeur. De plus, il arrive que la stratigraphie d’un site soit trop contrastée pour 
permettre de définir un module de sol équivalent. 


Dans notre étude, on va s’intéresser aux modules de déformation dynamiques et plus 
particulièrement au module de cisaillement dynamique G qui s'avère un paramètre 
essentiel pour tous les problèmes de déformation en géotechnique et les études 
particulières de comportement aux vibrations (dues aux tremblements de terre, au trafic, 
aux machines, aux explosions... ). 


Ill. Le module de cisaillement aux très faibles déformations, Go: 

Le module de cisaillement G est constant et maximal pour des déformations 
inférieures à la limite élastique (y < y£ = 107”). Il est donc noté Go où Gmax. En mettant 
au point des mesures de déformations très précises. Il a été possible de montrer que les 
modules de cisaillement dynamique et statique aux petites déformations étaient les 
mêmes. Par contre, pour des déformations plus importantes, le module dynamique est 
supérieur au module statique. En effet, d’après certains essais, la vitesse de chargement a 
peu d’influence sur le module de cisaillement pour des petites déformations.Il existe 
plusieurs méthodes pour la détermination du module de cisaillement qu’on peut les 
classer en catégorie : Méthodes statiques et méthodes dynamiques. 
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1. Méthodes statiques : 


Par essai statique, la vitesse de chargement (monotonique ou cyclique) est 
suffisamment lente pour que les effets de l’inertie soient négligés. Seuls les essais de 
laboratoire utilisent des méthodes statiques pour la détermination du module de 
cisaillement aux petites déformations. Ils sont devenus plus utilisés ces dernières années 
grâce à l'amélioration des méthodes de mesures de plus en plus fines des déformations et 
des forces. 


Notons toutefois que sur le terrain, les essais de plaques ou les essais pressiométrique 
permettent d'obtenir le module d'Young, mais les déformations mesurables sont 
supérieures à la limite élastique et donc ne correspondent pas au domaine de 
déformations considéré dans cette étude 


1.1. Essai triaxial 


L'essai triaxial permet d'obtenir le module de déformation longitudinale ou module 
d'Young noté E dans un domaine de déformations assez étendu, comportant notamment 
les déformations inférieures à la limites élastique (yê = 107”). Il est alors possible 
d'obtenir le module de cisaillement G grâce à la relation G = f(E) valable dans le domaine 
élastique: 


avec : v étant le coefficient de Poisson. 





0 >(1+v) 


capteur de force 


mesure de u (capteur) 
ou 


Contrepression mesure de Av (volumetre) 


eventuelle 
Figure 9.5 : Principe de l'appareil triaxial de revolution 





Figure 22: Principe de l'appareil triaxial de révolution 
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Remarque : 


La difficulté majeure des essais statiques ou monotones réside dans la mesure très 
fine des petites déformations. Effectivement, les erreurs apportées par exemple par 
l'appareillage utilisé deviennent d'autant plus significatives que les valeurs à obtenir sont 
faibles. Sur un essai triaxial par exemple, la friction du piston et les hétérogénéités aux 
extrémités des échantillons (contacts échantillon-base et échantillon-tête de chargement) 


participent activement à l'erreur respective des mesures des contraintes et des 
déformations. 


1.2.Appareil de torsion : 


L’essai de torsion permet d’évaluer directement les caractéristiques de 
cisaillement du sol étudié : t (contrainte de cisaillement) et y (déformation de 
cisaillement ou distorsion). Alors que dans les essais triaxiaux, les contraintes et les 
déformations appliquées sont longitudinales. 


Cet essai consiste à appliquer un moment de torsion T à un échantillon 
cylindrique et à mesurer les déformations engendrées par cet effort, l’étude du 
comportement du sol en contrainte-déformation est donc directe, comme pour l’essai 
triaxial. Malgré un appareillage et une procédure un peu complexe ainsi qu’une 
distribution non uniforme des contraintes et des déformations dans l’échantillon, cet essai 
présente les avantages suivants : bonne simulation du cisaillement simple ; accès aux 


petites déformations par le biais d’un matériel de mesure adapté et d’une bonne liaison 


entre les extrémités de l’échantillon, le socle de l’appareil (partie fixe) et la tête de 


€ ı Excitation de torsion 
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Tête rigide 


l’appareil (partie mobile). 
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Figure 23: schéma de l'essai de torsion 
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2. Méthodes dynamiques (sismiques) 


Les méthodes dynamiques sont associées à la propagation des ondes de 
cisaillement dans les sols qui permet d’obtenir le module de cisaillement dans la mesure 
où l’étude se limite aux petites déformations : la limite élastique n’étant pas dépassée, le 
comportement du sol est supposé élastique. On peut alors obtenir le module de 
cisaillement Go (ou Gmax) à partir de la vitesse de propagation des ondes de cisaillement 


Vs dans le sol considéré par la relation suivante : 


Go = PXV* , 
p : la masse volumique du sol exprimée en kg/m’ 
V; : la vitesse de propagation des ondes de cisaillement exprimée en m/s. 
Remarque : 


Cette relation du module de cisaillement en fonction de la vitesse d’onde, on va la démontrer 
mathématiquement par la suite. 


2.1.Méthodes sismiques de terrain 


Sur le terrain, les méthodes sismiques ou géophysiques, basées sur la propagation des 
ondes, sont très répandues. L’émission et l'arrivée du signal sont visualisées sur un 
système d’acquisition, ce qui permet d’obtenir la vitesse de propagation des ondes et d'en 
déduire le module de cisaillement ( Go = p X V.° ) dans la mesure où les déformations 
engendrées sont inférieures à la limite élastique. 


Globalement, on distingue les méthodes: où les instruments de mesure sont placés dans la 
masse du sol des méthodes où tout l*appareillage est en surface. 


2.1.1. Essais dans la masse du sol 


Le(s) récepteur(s) et le générateur peuvent être placés a des profondeurs donnée: soit en 
effectuant des forages préalables à l’intérieur desquels seront placées les cellules de 
mesure, soit, comme il l’a été récemment développé, en introduisant les cellules de 
mesure: dans des appareils de pénétration tels que le cône de pénétration CPT. 


L’échelle de mesures de 5 a 20 m est celle des structures de génie civil. Ces mesures 
peuvent être faites assez exceptionnellement à partir de la surface du sol ou, c’est le cas 
général, à partir des sondages ; on réalise l’ensemble des essais classiques sous les noms 
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de cross-hole, down-hole et up-hole. 


Enregistreur : Enregistreur ~ Enregistreur 
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Figure 24: Essais sismiques pour la caractérisation du sol 
a)Cross-hole b) Down-hole c) Up-hole 


2.1.1.1.L’essai cross-hole : 
Les essais cross-hole se font à partir de deux ou trois forages : 


v un forage dit forage d'émission ; 
v un ou deux forages de réception. 


La mesure de la vitesse de propagation de l’onde se fait entre deux forages parallèles 
distants de quelques mètres. L’impulsion est donnée dans le forage d’émission et est 
reçue dans le forage d’écoute à la même profondeur. Cette méthode permet de donner en 
fonction de la profondeur la valeur moyenne en plan des vitesses de propagation des 
ondes. 


L’essai cross-hole est très séduisant du fait de sa simplicité relative et de la qualité 
d'informations qu’il permet d’obtenir 


2.1.1.2.L’essai down-hole : 


La mesure de la vitesse de propagation de l’onde se fait le long d’un forage. L’émission 
du signal a lieu à la surface et la réception se fait à l’aide de capteurs placés dans le 
forage. Il est possible avec des sources d’énergies adaptées (frappe horizontale d’un 
massif posé en surface, par exemple) de donner à une forte proportion d’ondes de 
cisaillement. 


Les valeurs obtenues dans un essai down-hole correspondent aux caractéristiques du 
terrain au voisinage du forage, pour une direction verticale de propagation des ondes. 


H.Z/A.M — Département Génie Civil — Option BPC 50 


Ecole Mohammadia d'Ingénieurs- Projet de Fin d'Etude 2016-2017. 


Cette direction de propagation corresponde à celle couramment retenue en génie 
sismique. Théoriquement, avec un espacement resserré des capteurs, 1l est possible de 
détecter des couches de plus faible caractéristiques, même si celles-ci sont incluses entre 
deux couches plus résistantes. C’est un des avantages majeurs de la méthode. 


Dans le cas du dispositif down-hole , le signal est émis en surface ou dans un forage et 
réceptionné à divers niveaux dans un forage situé à proximité de façon que la plupart des 
trajets soient subverticaux. 


2.1.1.3.L”essa1 up-hole : 


L’essai est en tout point analogue à l’essai down-hole à la différence pres que l’émission 
est provoquée en fond de forage et la réception se fait sur des capteurs placés dans le 
forage, au-dessus du point d'émission. Par rapport à l’essai down-hole, l’exécution est 
plus délicate en raison de l’utilisation des dispositifs de frappe spéciaux et donc 
l’émission de l’onde de cisaillement est difficile. 


Généralement, et pour des raisons pratiques, on se limite à un ou deux points de réception 
par couche pour une même émission. 


Les essais down-hole et up-hole se font à partir de deux forages, ou à partir de la surface 
avec la réception dans un forage. Ces forages sont situés le plus près possible l’un de 
l’autre. On choisit un point dans un forage, de façon à avoir des trajets les plus verticaux 
possibles. 


Le calcul des vitesses nécessite une connaissance précise des distances, qui est obtenue 
par des mesures d’inclinomètres des forages. La précision qu’on peut attendre des 
mesures de temps, de distance et de masse volumique conduit à des modules déterminés 
en moyenne à 10 % près. Les écarts augmentent dans les terrains sans cohésion dans 
lesquels l’énergie des ondes de cisaillement semble s’amortir très vite. On est amené à 
choisir des distances émission- réception faibles, ce qui nuit à la précision relative des 
mesures. 


2.1.2. Essais en surface 


Dans ce cas, émetteur et récepteur(s) sont placés à la surface. Le problème majeur de ces 
techniques est qu’on ne maîtrise pas exactement les profondeurs pour lesquelles ces 
essais sont représentatifs ; est-ce une mesure ponctuelle ou une moyenne ? 


Les ondes de Rayleigh, encore appelées ondes de surface, se propagent dans des couches 
de sol dont l’épaisseur est égale à la longueur d’onde émise. La vitesse de propagation 
des ondes de Rayleigh est reliée à la vitesse de propagation des ondes de cisaillement Vs 
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par le coefficient de Poisson v. Toutefois, la relation simple Vg = 0,9 x V donne de bons 
résultats pour une gamme importante du coefficient de Poisson v (Richart et al, 1970). 


Emetteur de vibrations 


Systeme d'acquisition 
et analyse de spectre 


Oscillateur et Géophones 


Amplificateur 





Figure 25: Essai des mesures des ondes de surface (Rayleigh) 


La méthode communément employée pour l’analyse des ondes de Rayleigh est 
indiquée sur la figure au dessus. Les ondes sismiques sont émises par le biais de 
vibrations de fréquence fixée et sont reçues par deux géophones placés à une distance d 
l’un de l’autre. 


Le temps de propagation de l’onde entre les deux géophones est donné par : 


SE 


Tf) = F 


Où f est la fréquence de l’onde émise ; 
q est le déphasage entre les deux signaux reçus. 


Il est alors possible d’obtenir : 


La vitesse de propagation des ondes de Rayleigh ( VR(f) = ). d étant la distance 


T a ) 
entre les deux géophones), permettant ensuite d’accèder au module de cisaillement Go. 


La longueur d’onde Lr (Lr(f) = RD 





) permettant de déterminer la profondeur de la 


couche concernée par les ondes Pa 


Le profil des modules de cisaillement Go est donné par des essais à différentes fréquences 
qui représentent différentes profondeurs. 


Remarque : 


D’autres méthodes (analyse spectrale des ondes de surface (SASW) , essai de réfraction 
de surface , essai de réflexion de surface) reposent sur les mêmes principes. 
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2.2.Méthodes sismiques de laboratoire : 


Les essais de laboratoire peuvent être réalisés aussi bien sur des échantillons intacts que 
sur des échantillons reconstitués. 


La colonne résonante 


Les essais à la colonne résonante sont des essais reproduisant des séismes ; ils permettent 
de déterminer les caractéristiques du sol grâce aux sollicitations dynamiques sous petites 
et moyennes déformations. Par exemple, les modules sont déterminés à partir de la 
vitesse de propagation des ondes dans le sol. 


Ces essais consistent à atteindre la fréquence de résonance de l’échantillon en 
imposant des sollicitations cycliques longitudinales ou en torsion. Ils permettent de 
déterminer un lien entre les fréquences de résonance et les modules de déformation du 
matériau en petites déformations. 


Un échantillon cylindrique de sol est sollicité en vibrations à amplitudes faibles (flexion 
et torsion). Le sol est supposé élastique linéaire et l’échantillon est assimilé à une 
poutre finie élastique linéaire soumise à une excitation périodique à l’une de ses 
extrémités. Cette hypothèse est utile pour l’identification de sa fréquence propre 
correspondant au premier mode de vibration. 


Un générateur de fonctions transmet un signal à des bobines montées sur l’échantillon qui 
entraînent en déplacement les noyaux magnétiques. La fréquence du signal, mesurée 
précisément par un fréquencemétre, est augmentée progressivement pour atteindre la 
résonance, repérée, par exemple, à partir du maximum de la courbe fréquence 
accélération (l’accélération est mesurée par un accéléromètre) 
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Figure 26: Modes de sollicitation d'un essai à la colonne résonante 
L’échantillon est fixé à une extrémité tandis que l’autre extrémité est soumise à une 


excitation périodique (torsion ou flexion). Les vibrations se répétent à différentes 
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fréquences jusqu’à ce que la résonance soit atteinte. Les grandeurs mécaniques 
caractéristiques du sol (module d’Young, module de cisaillement, vitesse de propagation) 
sont déterminées à partir de la fréquence propre du premier mode de vibration. 


[Module] = [Masse volumique] x [Vitesse de propagation] * 


Or, par conception, l’éprouvette vibre avec la plaque motrice. La détermination des 
vitesses est de ce fait complexe. Il faut en effet réaliser des étalonnages pour déterminer 
les paramètres propres de l’appareil. Les formules suivantes constituent une procédure 
de la détermination de la vitesse de propagation de l’onde du cisaillement : 


Vs est obtenue avec la relation théorique 


Avec: 


f: fréquence de résonance, 


p : coefficient, qui dépend du moment d’inertie de l’éprouvette et du système de 


mise en résonance, déterminé par : = = p tan(b ) 
0 


mD? 
Avec: I = > 





e lo: moment d'inertie massique d’étalonnage, dépendant du système. Il se calcule 
par étalonnage à l’aide de l’éprouvette en aluminium de différentes inerties et de 
masses additionnelles à positionner sur le chapeau. Il est égal a 0,0039844 pour le 
système d’un système utilisé a LCPC, 

e I: Moment d'inertie massique de l’échantillon 


L : longueur de l’échantillon 
e D: diamètre de l’éprouvette. 
m : masse de l’échantillon. 


2.3.Comparaison laboratoire-terrain: 


Notons dans un premier temps que le prélèvement des échantillons in situ, malgré toutes 
les précautions prises durant cette étape, ne peut éviter un remaniement. La reconstitution 
des conditions en place est de la même manière une phase délicate des essais de 
laboratoire. 
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Les recherches visant à comparer les modules de cisaillement Go ou les vitesses de 
propagations des ondes de cisaillement Vs, mesurés par des essais en laboratoire et des 
essais in situ sont assez contradictoires. Il semble qu'il y ait toutefois une tendance à 
trouver des valeurs de module Go avec les méthodes de laboratoire plus faibles que ceux 
issus de méthodes de terrain. 


Un certain nombre d’études ont effectivement abouti à des écarts importants entre le 
terrain et le laboratoire. Stokoe et Richart (1973) ont comparé ces deux types d’essais 
avec la colonne résonante et des essais « crosshole ». Les vitesses de propagation des 
ondes de cisaillement étaient plus faibles en laboratoire (de 95% pour un sable silteux à 
70% pour un silt argileux ). 


Conclusion : 


On peut dire que la détermination des modules dynamiques semblent une tâche très 
délicate vu la difficulté de l’évaluation des caractéristiques des sols 


3. Détermination de module de cisaillement en fonction des vitesses 
d’ondes : 
On suppose que le milieu semi-infini, dont l’onde sismique se propage, est 


homogène et isotrope. On considère tout d’abord un élément infiniment petit de 
dimensions dx, dy, et dz. 


On considère également la variation des contraintes au niveau de chaque face de cet 
élément. Ces contraintes sont représentées par un ensemble des vecteurs orthogonaux. 





Figure 27: Elément volumique infiniment petit 
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L’équilibre de translation de cet élément est exprimé par la somme des forces appliquées 
parallèlement au chaque axe 


Par exemple, suivant la direction Ox , l’équation d’équilibre s*écrit sous la forme : 


(a, + ox dx) dydz — 0,dydz + es + Za dy) dxdz — t,,dxdz + Ca + 


X 
ð 
2 dz) dxdy — t,gdxdy (1) 
Des équations similaires à celle-ci sont exprimées suivant les directions Oy et Oz. 


En négligeant le poids propre de cet élément et en appliquant le 2°" loi de Newton on 
obtient suivant la direction Ox : 





00% , OTx ÔTxz a 
(= +72 + az) dxdydz = p(dxdydz) + (2) 
Donc : 


cd OT xy OTxz (3.1) 
Ot 0x dy OZ 





p 


De même pour les autres directions Oy et Oz : 


Ov 00, ÔtTyx ÔTy 
A 
Tear a 


ðw 00, OTzx  0OTzy 
— = 2 3.3 
ðt? dx dl dy Oz a 








p 


Avec u,v et w sont les déplacements suivants les directions Ox , Oy et Oz respectivement 


On peut exprimer le côté droite des équations au dessus en fonction des déplacements par 
les relations suivantes qui sont valables pour un milieu élastique : 


Ox = ÀE + 2GE, Try = Tyx = GYxy (4.1) 
Oy = ÀE +2Ge, Tyz = Tzy = GYyz (4.2) 
O, = le + 2Ge, Tri = Te = UY (43) 
G = 204v) (5.1) 
1 Ed (5.2) 


Qta- 2v) 
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Avec : v: coefficient de poisson 
A:Constante de Lamé 
G : module de cisaillement 


E : Module d YOUNG 


ml 


= Ex t Ey t Ez 


2016-2017. 


De plus, les déformations et les rotations peuvent êtres exprimés en fonction des 
déplacements ( Timoshenko and Goodier , 1951 ; Kolsky , 1963 ) : 


Y” Les déformations : 


E du dv A du i 
E dv B ðw A dv > 
E ðw B du A ðw e 
Ez ss 0z Vzx nn 0z Ox ( x ) 
Y” Les rotations: 
=S ðw ðv == du ðw == Ov ou 
20, = öy = Oz (7.1) 20) = à = En (7.2) 20, = Jx = dy (7.3) 


Avec 0, ,@y et wz sont les rotations autour de Ox , Oy et Oz respectivement . 


D’après les équations (3), (4) et (6) on obtient pour chaque direction les équations 


suivantes : 


Y 


% Suivant OX : 


2 


de 


— — == 2 
poa U+G)_+6Vu 


Y 


% Suivant Oy: 


2 de 
— = (4+ + GVév 
Poe dy 


UV 


% Suivant OZ : 


2 


ð w de a 
pe A0 + GV°w 


Oz 


H.Z/A.M — Département Génie Civil — Option BPC 


(8.1) 


(8.2) 


(8.3) 


57 


Ecole Mohammadia d'Ingénieurs- Projet de Fin d'Etude 2016-2017. 


; e A: 07 
Avec V° : est Laplacienne dans les coordonnées cartésiennes V? = Far taz 
Les équations (8.1) , (8.2) et (8.3) sont les équations de mouvement d’un milieu 
élastique , semi-infini , homogène et isotrope 


On différencie les équations (8.1), (8.2) et (8.3) respectivement par rapport x , y et z 
puis on effectue la sommation on obtient donc : 


2 


0% E 
ae ARN E 


Pa 
Orona: a = v VE (9) 
? : ’ 2 (2+2G) 
D'ou la formule de la vitesse d’onde : v,* = eo 


ve: c’est la vitesse d’onde de compression dans un milieu continu 


L’autre solution des équations du mouvement sont obtenus en différenciant (8.2) par 
rapport à z et l’équation (8.3) par rapport à y. Puis , on élimine € en effectuant la 
différence des 2 équations , on obtient donc : 


= Gv? (2 m T 
E Oy ðz SA) 


On utilise par la suite l’expression de la rotation w,équation (7.1) : 


ð? fs 2 
Par Nay oz 








0 w. o 
P o = GV ©; 
Orona: 
0 w. en 
TA = v V O, (11) 


Les mêmes expressions peuvent être obtenues pour wy et w, ce qui implique : 


Us = D 


On obtient par la suite la fameuse équation du module de cisaillement en fonction de la 
vitesse d’onde de cisaillement : 


G = px? 
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Chapitre V : interaction sol-structure (ISS) 


I. Généralités : 


La réponse de la structure à un mouvement sismique dépend de la nature de ce dernier, 
des propriétés dynamiques de l’ouvrage et de celles du sol. Il y a alors une interaction 
entre la structure et le milieu extérieur. 


Donc, on peut dire que L’interaction sol-structure est une discipline de la mécanique 
appliquée s’intéressant à l’analyse des structures soumises à des charges dynamiques en 
tenant compte du comportement du sol de la fondation. Les effets de l’interaction sol- 
structure (ISS) sur la réponse sismique n’ont été sérieusement pris en considération 
qu'après le tremblement de terre de 1971 à San Fernando et au début de la construction 
nucléaire en Californie. L’analyse sismique d’une structure tenant compte des propriétés 
locales du site, diffère de celle considérée encastrée à sa base De ce fait, 1l est 
particulièrement important de considérer l’ISS dans les zones sismiques où la réponse 
dynamique des sols peut changer la réponse des structures soumises à l'excitation 
sismique. De plus, les sites très meubles (sol mou) recevant des structures rigides et 
massives peuvent changer les caractéristiques dynamiques des ces derniéres de maniére 
significative. 


Et pourtant, le règlement parasismique marocaine RPS 2011 tient compte les effets de 
PISS d’une manière implicite en proposant les catégories de sites en fonction des 


caractéristiques du sol qui les constituent, chaque catégorie de site est associée à un 
spectre de réponse élastique obtenu à partir des propriétés du site et de la structure. 


Remarque : 


Les raisons pour lesquelles l’interaction sol structure n’est pas plus souvent prise en 
considération dans les études sont essentiellement liées à la méthode et l’outil, ainsi que 
le manque de données de sol. De plus, la prise en compte de l’interaction sol-structure 
complexifie beaucoup la modélisation de la structure. C’est pourquoi d’une façon 
conservatrice, l’interaction sol-structure est souvent négligée lors des analyses sismiques. 


II. Influence de la nature des sols sur le mouvement sismique 
L’observation montre que l’intensité avec laquelle un séisme est ressenti en un lieu 


donné dépend dans une large mesure de la nature des terrains traversés par les ondes 
sismiques et des conditions locales. 
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Il est constaté souvent que les ouvrages édifient sur un sol mou subissent des 
dommages plus importants que ceux situés sur un sol rocheux. Ce phénoméne peut 
expliquer par le fait que les couches de sol très meubles se comportent 
comme un oscillateur qui amplifie l’excitation appliquée à la base par le rocher. Il 
convient donc d’adapter le spectre de réponse à la nature du sol. Les sols meubles filtrent 
les fréquences supérieures de séismes et ne laissent passer que les fréquences basses qui 
sont généralement voisines de celles des constructions. 


Les spectres de réponse moyens sur la figure au dessous calculés à partir 
d’enregistrements réels obtenus sur des sols rocheux, argileux, raides, de profondeur 
moyenne ou des dépôts alluviaux profonds, montrent la nécessité qu'il y à étudier le 
comportement des sols sous chargement dynamique pour parvenir, en zones sismiques, à 
une conception fiable des ouvrages qui y sont fondés, c'est en particulier le cas pour les 
ouvrages importants comme les centrales nucléaire, les barrages, les tours, ...etc 


Accélération epectrele 
3 Accélération maximale du s 


A 


Argile de consistance molle À moyenne et sable 
| Sals prolonds sans cohésion 


| Sols dé consislance raida 





Figure 28: Allures des spectres d'accélération moyenne pour différents sols 


II. Nécessité de modéliser le sol 


Lors d’un tremblement de terre, les ondes sismiques se propageant dans le sol mettant en 
mouvement les fondations des ouvrages en les secouant principalement horizontalement. 
Accéléré a sa base, chaque bâtiment est soumis à des forces d’inertie auxquelles sa 
structure doit résister. Pour maintenir son équilibre, l’ouvrage exerce sur le sol 
environnant des efforts importants. Si ce sol est très résistant, le mouvement initial de 
l’appui n’est pas modifié. Il n’y a pas d’interaction et pour assurer la tenue de l’ouvrage, 


le concepteur peut restreindre son analyse à la seule structure ; Si par contre le sol est de 
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résistance plus faible, le support va se déformer et le mouvement de la fondation s’en 
trouve modifié. Nous avons affaire à un phénomène d’interaction sol-structure et le 
concepteur doit alors inclure le sol et parfois les ouvrages environnants dans son analyse. 


IV. Limitation de PISS : 


Tout d’abord, la modélisation la plus fréquente de PISS suppose que la structure est 
fondée en surface sur un sol homogène (solide élastique) et horizontal, avec des 
propriétés mécaniques connues et constantes au cours du séisme, et que le mouvement 
sollicitant la structure soit des ondes de volume se propageant verticalement à partir d’un 
substratum horizontal. 

En réalité, le sol est hétérogène et les caractéristiques dynamiques varient au cours du 
séisme, la structure peut être aussi sollicitée par des ondes de surface, et la fondation, 
voire une partie de l’ouvrage, est le plus souvent disposée à une certaine profondeur. 
Enfin, le solide élastique est supposé avoir le même module en traction qu’en 
compression, ce qui se concilie mal avec un sol ; ce module est supposé le même en tout 
point du sol, ce qui est en contradiction avec la réalité. 

L’ensemble de ces incertitudes fait que le calcul de PISS reste encore approximatif et 
donc une certaine prudence s’impose à l’exploitation des résultats obtenus. Par ailleurs, 
ces incertitudes conduisent à paramétrer les caractéristiques du sol (voir la partie de la 
dynamique du sol ) et à prendre l’enveloppe des résultats. 

Dans telles conditions, l’ingénieur est souvent amené à raisonner de manière différente 
plus intuitive et pragmatique, pour tenter de trouver une solution satisfaisante au 
problème si complexe de l'ISS. 

On est ainsi amené à concevoir des modèles dans lesquels le sol est représenté par des 
ressorts ou des éléments finis. 





Figure 29: Types de modèles 
a) Encastrement parfait b) Ressorts c) Eléments finis 
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V. Les conséquences de l’interaction sol-structure 
Le problème de l’interaction sol-structure peut être décomposé en deux 
ordres classiquement dénommés : interaction cinématique et interaction inertielle. 


1. L’interaction cinématique: 


Cette interaction traduit l’incompatibilité entre le champ d’onde incident et les 
mouvements possibles de la fondation qui est considérée sans masse, vu la discontinuité 
des rigidités (l’excavation et l’inclusion de la fondation ont changé les caractéristiques). 
Ce qui empêche la fondation de suivre les mouvements imposés par le sol. le champ 
d'ondes incident est réfléchi et diffracté par la fondation et donc modifie le mouvement 
total du sol au voisinage de celle-ci et réduire l’énergie pénétrant dans la structure. Donc, 
cette interaction cinématique induit un filtrage parfois significatif du mouvement 
sismique lors de son transfert au bâtiment mais occasionne des rotations d’ensemble de la 
fondation dont les effets sur les ouvrages doivent être étudiés. La réponse de la fondation 
sans superstructure au mouvement sismique fournit au concepteur les efforts d’inertie 
nécessaire au dimensionnement de l’ouvrage. 


Remarque : 


Dans le cas de fondations superficielles et d’ondes à incidence verticale, il n’y a pas 
d'interaction cinématique et le mouvement de la fondation sans superstructure est égal au 
mouvement du sol sans ouvrage. 





Figure 30: Interaction cinématique entre une fondation superficielle et une onde plane 
en incidence oblique. 


2. L’interaction inertielle: 


Cette interaction réside dans le fait d’inclure la masse de la structure, qui va induire des 
efforts d’inerties et des moments supplémentaires à la base de la structure et donc dans le 
sol, ce qui aura pour conséquence de modifier encore une fois le champ de déplacement. 
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Par ailleurs, le mouvement induit sur la fondation qui développe des oscillations de la 
superstructure et donc donne naissance à des forces d'inertie qui sont retransmises à la 
fondation sous forme de forces et de moments. 


Remarque : 

Le dimensionnement de la fondation doit tenir compte de ces deux composantes de 
l'interaction. (Les effets de l’interaction cinématique sont plus difficiles à évaluer 
rigoureusement que les effets inertiels.) 


VI. Modélisation du sol en tenant compte de l’interaction sol 
structure 


Il existe deux méthodes classiques pour tenir compte de l’interaction sol-structure dans 
une analyse dynamique : 
% La méthode des éléments finis. 


Y 


% La méthode de la sous-structure. 


1. Modélisation par des ressorts : 


La manière la plus simple de prendre le sol en compte est le représenter par des ressorts 
reliant un ou plusieurs nœuds à une base rigide, à laquelle on impose le mouvement. 
Dans le cas d’un modèle plan, une semelle sous point d’appui isolé est représentée par 
deux ressorts agissant à la translation et un ressort à la rotation. Sous un radier, le sol est 
modélisé par un ressort horizontal et un ressort vertical en chaque nœud 


al Kn 


EJ 
Oa 
TOT 
IN 
RA 
de 
VS 


WIN 


K, 





VI LI III IIS, PILITI III, 


Figure 31: Modèle du système barres-ressorts 
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Remarque : 


Il y a pas mal des méthodes qui calculent les rigidités et également les amortissements 
pour chaque ressort, qu’on va les expliciter par la suite. 


2. Méthode des éléments finis : 


La modélisation par éléments finis est la plus directe mais souvent aussi la plus 
complexe, en particulier à cause de son temps de calcul. Le principe de cette modélisation 
est de représenter le sol par les couches successives d’éléments empilés sur une base qui 
représente le soubassement rocheux. Chaque couche est caractérisée par un module 
d'élasticité, un coefficient de poisson, une masse volumique et un coefficient 
d'amortissement. Ainsi elle nécessite l’utilisation des accélérogrammes en champ libre et 
à la base de la structure 


Eléments barres 


Couches successives 


Soubassement 
rocheux 





Figure 32: Modélisation par éléments finis 


Le mouvement est imposé a la base et on fait l’hypothèse qu'il se propage dans le sol 
sous forme d’ondes de cisaillement se déplaçant uniquement à la verticale. 
Le mouvement de la base (située en profondeur) est différent de celui de la surface du 
sol (en champ libre).Il n’est pas connu à priori et par conséquent on doit étudier au 
préalable le sol en absence de toute construction pour déterminer le mouvement de la 
base qui produit en surface le mouvement du champ libre défini par l’accélérogramme 


donné. 


Œ Suite aux difficultés que présente cette méthode et au manque des données on a 
opté pour la méthode de la sous-structure. 
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VII. Les amortissements : 


1. Les mécanismes d’amortissement : 


L’amortissement est une propriété dynamique d’un système mécanique qui s’oppose aux 
déformations et aux déplacements. 


Lors de l’interaction sol-structure, une partie d’énergie communiquée à la construction 
est renvoyée dans le sol et dissipée par deux mécanismes distincts : 


> Atténuation par radiation, appelé également amortissement géométrique ou 
pseudo-amortissement, car le phénomène est purement géométrique ; les ondes 
renvoyées ayant un front sphérique, mettant en vibration un volume de terrain de 
plus en plus grand. 

> Amortissement interne du sol appelé également amortissement matériel qui 
traduit la dissipation en chaleurs lors de ses déformations qui est due à l'élasticité 
imparfaite du sol. 


Ces deux types d’amortissements se produisent lors de la propagation de toute énergie 
mécanique dans le sol. 


Remarque : 


L”amortissement est appelé géométrique, car il ne dépend que de la forme de la 
fondation, des propriétés du sol et la fréquence du mouvement. 


2. Caractérisation de l’amortissement : 


L’amortissement du système sol-fondation est obtenu en employant les formules 
classiques qui sont indépendantes de la fréquence. On calcule les rigidités du sol 
indépendantes de la fréquence pour une première approximation par la méthode de 
Newmark-Rosenblueth ,d’une part, on peut calculer les rigidités et l’amortissement 
géométrique d’une manière précise, d’autre part, dans le cas des fondations circulaires ou 
assimilées sur un milieu semi infini, selon la méthode de Deleuze. 


Cependant, pour notre étude, l'analyse a été effectuée en considérant un comportement 
linéaire du sol avec un amortissement visqueux, ce qui ne sera pas applicable dans le cas 
de forte sismicité où les non linéarités de sol peuvent jouer un rôle décisif en modifiant 
l'état du sol autour de la fondation et en provoquant un fort amortissement de l'énergie 
sismique injectée dans la structure. 


Dans le cas de sollicitations sismiques faibles voire modérées, un modèle élastique avec 
un amortissement visqueux peut être représentatif du comportement du sol. 
L’amortissement est pris en considération par la formulation de Rayleigh 
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L”influence de la plasticité du sol sur l’interaction sol-structure a montré que la plasticité 
du sol induit deux effets, à savoir un amortissement supplémentaire dû à la dissipation 
par déformations plastiques et une réduction des « fréquences propres » du système sol- 
fondation due à la réduction de la « rigidité » induite par la plasticité. L’ampleur de 
l’influence de la plasticité dépend de son extension dans le massif du sol, qui dépend de 
l’amplitude du chargement, de son contenu fréquentiel et des fréquences propres du 
système sol-fondation structure. 


Et comme on a supposé que le milieu est semi-infini (demi-espace) élastique homogène, 
donc on aura une surévaluation des valeurs de l’amortissement radiatif, ce qui fait que 
l’amortissement de chaque ressort sera réduit et pris égale à la somme : 


amortissement du ressort = (amortissement radiatif ou 
géométrique)+(amortissement interne du sol) < 30% 


L”amortissement interne de frottement du sol est égal a 5 %, quels que soient le mode et 
le sol. 
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Chapitre VI : Les méthodes de calcul des raideurs et 
d’amortissement des ressorts 


Ce chapitre a pour but de résumer les différentes formules applicables pour déterminer 
les raideurs de fondations superficielles que l’on peut trouver dans les différents ouvrages 
abordant l’interaction sol-structure. 


Une première partie consiste à faire l’inventaire des méthodes les plus utilisées afin de 
déterminer les raideurs et amortissements des sols. La seconde partie contient un 
comparatif des différentes méthodes. 


I. Méthode de Newmark-Rosenblueth [1971] 


Cette méthode est indépendante de la fréquence de la structure et permet d’estimer et de 
contrôler les valeurs de raideurs et des amortissements du sol. Rosenblueth a introduit 
une masse virtuelle de sol, liée à la fondation du bâtiment. Cette masse ne représente pas 
une masse de sol « attachée» à la fondation mais constitue un artifice mathématique 
permettant de mieux prendre en compte de la réponse de la fondation. L’ajout de cette 
masse ne modifie pas radicalement la réponse dynamique de la structure. 


Les raideurs du sol ainsi que les amortissements sont donnés par le Tableau 4 : 


Facteur Raideurs K 
Hauteur du d'amortissement 2 
Mouvement ; i J Fondation . y 
prisme de sol | relatif ou critique À i Fondation rectangulaire 
circulaire 
n 
Vertical 0,27VA 


: 32(1 — v)Gr, 
G 
RA KE 
Rotation autour ebay 
de Paxe vertical 0,25VA = ——Gf, (a? + b2)VA 
(torsion) 4 


Tableau 7: les différents raideurs et facteurs d'amortissement selon la méthode de 
NEWMARK-ROSENBLUETH 





Avec : 
À : aire de la fondation 
a : dimension parallèle à la direction du séisme 
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b : dimension perpendiculaire à la direction du séisme 

Mb : masse du bâtiment 

Ms =A.H.p : masse de sol additionnelle ayant même aire que la fondation du bâtiment 
H : hauteur du prisme de sol 


2+H} Méx. | pon n I 
Jp = M; E + 7) inertie massique du bâtiment par rapport à la fondation 
— a? +Híg Hs : $ a Sri ` 
Js = Mo NT E ) inertie massique de la masse additionnelle de sol par rapport à 


la fondation 
ro : rayon d’une fondation circulaire 
v : coefficient de poisson 


Bx, Pz et Po sont données par les abaques suivants : 


ooo + tit | | ||| 
>= tp titi, | | | 
o ABU S ION | | |. 
A ESA TT 


_ 





T 


0.1 0.2 03 040.5 1.0 2 E 4 5 10 a/b 


Figure 33: Abaques donnant les coefficients de Newmark 6,, 6, et 6p 


Les valeurs des fy, z et Po Seront constants tout au long de notre étude, car les semelles 
utilisées pour nos cas seront des semelles carrées a=b alors : 


e 65,=0,97 
° bz = 2,2 
o Bo = 0,47 


La valeur retenue du facteur d’amortissement est la moitié de la valeur ci-dessus à 
laquelle on ajoute l’amortissement du sol de 5% : 


E=2+5% < 30% 
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ll. Méthode de Deleuze : [1967] 


Dans une analyse de l'interaction sol-structure, les coefficients dynamiques de la rigidité 
du sol à l'origine dépendent de la fréquence, selon la méthode de Deleuze, on peut 
calculer les rigidités et l’amortissement géométrique d’une manière précise. Cette 
méthode peut être appliquée dans le cas où les fondations, circulaires ou assimilées, 
peuvent être considérées comme superficielles et suffisamment rigides. Le sol est 
représenté par un système de ressorts qui relient le ou les nœuds représentatifs de la 
fondation de l’ouvrage au repère de sol en champ libre. 


Les coefficients de raideurs et les amortissements géométriques en fonction des 
coefficients de transmittance f sont donnés dans le tableau suivant : 


Mouvement Raideurs Amortissements géométriques 
fvi o= frz 


Vertical ky = Gro 55 Ny = 
2fyı 


få + fh 


Horizontal kn = Gro 5557 Ny = — 
DU HE | GG 2fu1 
= 3 fri Fo 

Balancement kg = GT — 3 No = — 
Fr + JR 2 fR1 


Tableau 8: les différents raideurs et facteurs d'amortissement selon la méthode 
DELEUZE 





Deleuze a établi des tableaux donnant des coefficients de transmittance fv (translation 
verticale), fx (translation horizontale) et fo (rotation autour d'un axe horizontal, 
balancement). Pour établir ces coefficients Deleuze a utilisé les hypothèses suivantes : 


+ Le sol est un solide semi-infini, élastique et homogène, isotrope et limité par un plan. 

+ La réaction de sol sous la fondation est linéaire. Les déplacements du nœud de la 
fondation sont assimilés a ceux du centre d'inertie de la fondation, la rotation de la 
fondation est celle du plan moyen de la surface d’appui 

+ L’amortissement géométrique correspond au rayonnement des ondes dans le sol 


Les tableaux suivants illustrent les différentes valeurs des transmittances ( vertical, 
horizontal et rotation) : 
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fv1 fv2 fv1 fv2 fv1 fv2 


O 026751 0 | 0,2026 | 0 | 0,1351 | 0 


aa a 


Tableau 9: Coefficients de transmittance verticale f, du sol 





fH1 fH2 fH1 fH2 fH1 fH2 fH1 fH2 


o [oz | o [oz o Joze) o ozs] o 


Tableau 10: Coefficients de transmittance horizontale f; du sol 


A a a a 





jaa] os | oa | vos 
ES fR1 fR2 fR1 fR2 fR1 fR2 fR1 fR2 
o [oss| o [o322| o [osa| o [ozs o 


Tableau 11: Coefficients de transmittance en rotation f, du sol 
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Remarque : 


e Lorsque les valeurs numériques des paramètres ne figurent pas dans les tableaux 
une interpolation linéaire est effectuée à partir des valeurs les plus proches. 

e Pour étudier des différentes cas, nous avons besoin de calculer à chaque fois les 
rigidités, et par conséquent, nous devons interpoler à chaque fois le coefficient de 
poisson v et le paramètre ao, ce qui semble une tâche délicate en raison de la 
diversité des cas et des rigidités, c’est pour cela, que nous avons programmé sur 
EXCEL à l’aide de VBA un programme qui permet d’interpoler afin d'optimiser 
nos calculs. 


Avec ro : rayon du cercle équivalent d’une fondation circulaire ; 1l est calculé comme suit : 


e Pour les raideurs en translation : l’aire du cercle équivalent est égale à laire de la 


i ba 
surface d'appul : rọ = ee 


e Pour les raideurs en rotations (balancement) : le moment d’inertie de la surface du 
cercle équivalent est égal au moment d’inertie de la surface d’appui de la fondation, 
rapporté a l’axe perpendiculaire a la composante passant par le centre de gravité de la 
surface, soit 


4 |baÿ 
v en balancement : rọ = pen 


4 |ba(a? +b?” 
v en torsion : y = | paata 


do = 2N f To £ : Paramètre utilisé dans les tableaux des coefficients de transmittance. 


En pratique cette méthode est applicable si ay < 2 . Celle-ci n’est applicable que si le 
rayon de fondation est suffisamment petit devant la longueur d’onde dans le sol. 


L”amortissement s’obtient avec les formules suivantes : 


1 

v Vertical : č, = >v + 5% < 30% 

l 1 
v Horizontal : ¿y = Mu + 5% < 30% 
l 1 

Y Rotation : čo = >e + 5% < 30% 
Ces formules donnent les raideurs et les amortissements réduits d’un ressort unique. Si le 
modèle ISS est celui des ressorts répartis, les raideurs correspondantes sont obtenues en 


privilégiant la direction horizontale et la direction verticale, dont les déplacements sont 
issus d’une rotation : 
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Y 


; : ; KH 
% direction horizontal = 


+ 


é : A Kọ 
% direction vertical = 


Avec: S : surface de la fondation et I : inertie de la fondation 


a : dimension parallèle à la direction du séisme ; 

b : dimension perpendiculaire à la direction du séisme ; 

Ro : rayon du cercle équivalent d’une fondation circulaire (voir tableau 11.2) ; 

p : masse volumique du sol ; 

G : module de cisaillement dynamique du sol ; 

f : fréquence du mode de vibration fondamental dans la direction de séisme considéré, 


IM. Guide du SETRA: 


Pour les semelles superficielles, dans le cas de structures simples à peu de degrés de 
liberté et des sols de stratigraphie régulière, l’interaction sol/structure peut être modélisée 
à l’aide de raideurs (ressorts K) et d’amortisseurs (amortisseurs C), pour une fondation 
circulaire équivalente reposant sur un demi-espace élastique, sont données ci-après : 


Fondation réelle Fondation circulaire équivalente 





Figure 34: Définition des paramètres utilisés dans les formules du guide SETRA 


Raideurs Amortissements 


_ 4GR = 
GR 


ko 


on 0,3 /Bo 
Co = ko 
31v) _ v) 1+B% 


Ie 0,3 ES R AUDE Bo = _ 3(1— n 
1+B9 Vs 8pR 


Tableau 12: Les a raideurs et amortissement selon le guide du SETRA 


7, Avec By = E 
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Vis-à-vis des raideurs associées aux degrés de liberté de translation, la fondation 
circulaire offre la même section que la fondation réelle soit mR? = a x b. Vis-à-vis de 
chaque degré de liberté de rotation, la fondation circulaire équivalente est celle qui 
présente la même inertie, 


| axb? : 
soit TR* = 2. Pour la rotation autour de l’axe Ox 
x 


bxa’ . 
et TR* = — pour la rotation autour de l’axe Oy. 


SCS 


e 


G le module de cisaillement 


% 


e 


v le coefficient de poisson 


% 


e 


p la masse volumique du sol 


SCS 


e 


Io et lo les moments d’inertie de l’ouvrage pour le balancement et la torsion 
respectivement 


+ 


Vs la vitesse des ondes de cisaillement 


IV. Comparaison des méthodes : 


Afin de choisir la méthode convenable à notre étude, on doit comparer tout d’abord ces 
méthodes tout en présentant les avantages et les inconvénients de chacune dans le 
tableau ci-dessous : 


Les méthodes 


NEWMARK- Calcul simple et rapide BOTO OUCIE | o 
ROSENBLUETH Fondation circulaire et Nécessite de connaitre la masse et l’inertie de la 
rectangulaire Structure 
Dépend de la fréquence Monocouche 
du mode fondamental de Calcul plus long car nécessite la détermination 
l’ouvrage de coefficient de transmittance. 
Nécessite le calcul d’un rayon équivalent dans le 
DELEUZE cas d’une fondation rectangulaire 


do = 2Tf To E doit être inférieur à 2 ce qui 


limite l’utilisation de cette méthode. 


Calcul simple et rapide Monocouche 


Guide du SETRA Nécessite le calcul d’un rayon équivalent dans le 
cas d’une fondation rectangulaire 





Tableau 13: Avantages et inconvénients des différentes formules de calcul de raideur 
et d'amortissement 
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On constate que ces méthodes ne sont applicables que pour des monocouches. Or dans la 
réalité il arrive très souvent que l’on ait affaire à un sol possédant plusieurs couches. Le 
PS 92 donne des formules pour calculer des caractéristiques équivalentes pour un profil 
stratifié dont les caractéristiques mécaniques varient peu d’une couche à l’autre. 





Pour la masse volumique du sol : p; = zor 
ad 2 GiHi 
Pour le module de cisaillement du sol : G; = SH; 


Avec Hi, pi et Gi les paramètres relatifs à la couche 1. 


D’après l’AFPS, les formules d’homogénéisation des couches de sol pour revenir à un 
modèle simple monocouche sont à utiliser avec précaution. En effet, ces méthodes sont 
limitées pour des profils de sols ou le contraste des propriétés au sein des différentes 
couches n’est pas important (rapport des Vs entre 2 couches compris entre 0,5 et 2). 
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Chapitre VII : Modélisation avec prise en compte de 
l’interaction sol-structure 


L Introduction : 


Pour la plus part des études sismiques, on considère la structure ancrée à sa base 
dans un milieu infiniment rigide, sans tenir en compte la souplesse et la flexibilité du 
sol et les éléments de fondation, et sa capacité à dissiper une partie de l’énergie sismique 
qui peut conduire à des résultats sensiblement différents. 


Ce travail consiste à calculer la réponse de la structure (avec et sans ISS) en terme 
des périodes fondamentales pour différents types de sites, ainsi que l’effort tranchant à la 
base de la structure; les déplacements et les déplacements inter-étages dans le plan du 
portique; puis on va étudier l’influence de la hauteur du bâtiment sur les paramètres de 
la réponse dynamique (période et déplacement). 


En un premier lieu, nous traiterons la modélisation d’un portique étagé (4 cas) 
sans tenir compte de l’effet d’interaction sol-structure dans les différents sites par 
l’utilisation de spectre de réponse associé à chaque site. 


En seconde lieu, la même structure sera modélisée en tenant compte de l'effet 
d'interaction sol-structure, respectivement pour les différents sites. 


Les réponses sismique avec et sans interaction sol-structure seront étudiées et 
comparées à l’aide du logiciel de calcul de structure ETABS, dans lequel les effets du sol 
sur les fondations, ont été représentés par des ressorts linéaires, cette analyse numérique 
sur le logiciel va nous permettre de déterminer les périodes pour chaque mode propre, les 
déplacements latéraux, et déterminer les spectres de réponse, l’effort tranchant à la base 
et les différentes déformées. 


Remarque : 


On va expliciter par la suite la manière dont on a modélisé nos structures sur le logiciel. 


IL Classification des sites selon le RPS 2011 : 


La classification des sites est déterminée sur la base des paramètres géotechniques 
résultant de la reconnaissance des sols qui le constituent. Selon RPS 2011, le sol est 
classé en 5 sites dont on va s’intéresser que pour les 4 premiers : le sol Sı, le sol 
S2, le sol S3, le sol S4. 
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Le sol Sı (site rocheux): roche ou autre formation géologique sain à moins de 
3 m de la fondation caractérisée par une vitesse moyenne d’onde de 
cisaillement Vs > 760 m/s. 


Le sol S2 (site ferme) : rocher altéré , sols cohérents très raides,sols 
pulvérulents très denses , marnes ou argiles très consolidés (la vitesse des 
ondes de cisaillement Vs atteint au moins 360m/s ) 

Le sol S3 (site meuble) : dépôt de gravier ou sable moyennement compacts ou 
argile moyennement raide dont les épaisseurs varient de quelques dizaines à 
plusieurs centaines de mètres, où la vitesse des ondes de cisaillement est d’au 
moins 180m/s 

Le sol Ss (site mou) : ce sol est constitué par des dépôts de sables lâches, sans 
ou avec présence de couche d’argile molle, Ou dépôt d’argile molle à 
moyennement raide, ou tout sol de profondeur de plus de 3 m ayant un teneur 
en eau w> 40% et indice de plasticité PI > 20. où Vs est inférieure à 180m/s 
sur les vingt premiers métres, ce sol peut aussi être constitué de sols cohérents 
de faible raideur. 


N.B : 


Le site ferme va être éliminé par la suite en raison que les résultats obtenus pour ce site 
seront proches à ceux du site rocheux, et par conséquent, notre étude sera restreinte sur 
trois sites (rocheux, meuble et mou). 


III. Présentation de la structure : 


Dans un premier temps, nous traitons un cas simple qui consiste en un portique en béton 
armé de 6,0 m d'ouverture. Les poteaux de 3m de hauteur, ont une section de 0,3x0,3 m2. 
La structure repose sur des semelles superficielles en béton de section carrée qu’on va 
les pré-dimensionner dépendamment des charges appliquées. Ensuite, on va généraliser le 
travail pour une large gamme de structures bidimensionnelles comportant plusieurs 
étages et plusieurs travées. Le sol est supposé homogène avec un comportement élastique 
linéaire. 


Ces structures ont été choisies pour effectuer une étude paramétrique pour évaluer le 
comportement dynamique d'une structure de bâtiment en vibrations forcées, 
Chaque structure est soumise aux charges dynamiques déduites à partir du spectre de 
réponse du code Marocain RPS2011 . 


La structure est classée dans les bâtiments de classe III et supposée être implantée en 
zone de moyenne sismicité (zone IT). 
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Dans le cas sans ISS, on utilisera l’hypothèse de l’encastrement parfait a sa base. Il s’agit 
donc du modèle utilisé classiquement pour les bâtiments courants. 
Les propriétés mécaniques de structure donnent comme suit ; module d’élasticité : 
Ep= 32000Mpa, coefficient de poisson v= 0.2. 


Les charges appliquées sur chaque travée : 


v Charge permanente G= 40 kN/m 
v Charge d’exploitation Q=10 kN/m 


La réponse de la structure est évaluée en considérant les spectres de calcul 
associés aux différentes catégories de site qu’on va les insérer sur le code du calcul 
ETABS afin d’entamer nos calculs. 


Caractéristiques des différentes catégories de site : 


Le module dynamique mesuré par les essais in-situ est une mesure sous sollicitations 
faibles, pour lesquelles les modules d’élasticité E et de cisaillement G sont élevés. On sait 
que sous sollicitations fortes (déformation induite par le séisme), le module de 
cisaillement G à prendre en compte est plus faible que celui sous sollicitations faibles. En 
raison, de tenir compte des déformations engendrées par les ondes sismiques. 


Dans la modélisation de PISS, on doit tenir compte la diminution du module de 
déformation avec l’augmentation du niveau de déformation. Pour cela, l’Eurocode 8-5 
propose de déterminer cette valeur à partir du module tangent Gmax, auquel on applique 


un coefficient réducteur en fonction de l’accélération du sol (tableau 11). 


Gs max 


Tableau 14: Coefficients de réduction moyen (+- écart type) pour le module de 
cisaillement G 





e Gmax est le module de cisaillement moyen à faibles déformations (mesuré par les 
essais in-situ) qu’on va le calculer d’après la relation : G = p Xx Vz 
e S: Paramètre caractéristique de la classe de sol défini dans l’EN 1998-1:2004 
Agr: est l'accélération au rocher 
Oo a a A per 


g yı: coefficient d'importance 


Le coefficient yı est donné dans le tableau suivant ($4.2.5 NF EN 1998-1 de sept 2005) : 
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Catégorie Catégorie de bátiment Y 
d' mportance 


Bâtiments d'importance mineure pour la sécurité des 
personnes, par ex. les bâtiments agricoles 


Bâtiments dont la résistance aux séismes est importante 
compte tenu des conséquences d'un effondrement, par ex. les 
écoles, les salles de réunion, les institutions culturelles etc... 
Bâtiments dont l'intégrité en cas de séisme est d'importance 
vitale pour la protection civile, par ex. les hôpitaux, les 
casernes de pompiers, les centrales électriques etc.... 
Tableau 15: Coefficient d'importance en fonction des catégories du bâtiments selon NF 
EN 1998 





Pour notre cas, on est dans la catégorie des batiments courants (catégorie IT) donc 

y =1 
Par ailleurs, 1l est nécessaire de connaître la classe du sol (décrivant les profils 
stratigraphiques et prenant en compte l'influence des conditions locales de sol sur l'action 


sismique). On définit ainsi les classes de sol A, B, C, D et E (cf tableau 3.1 NF EN 1998- 
l de septembre 2005). 


Classe de sol 
l 





Tableau 16: le paramètre caractérisant la classe de sol type 1 (forte sismicité - zone 5) 
selon l'EC8 





Tableau 17: le paramètre caractérisant la classe de sol type 2 (faible sismicité - zones 1 
à 4) selon l'EC8 


Pour notre étude, on est dans une zone de faible sismicité donc on va travailler avec le 
2 tableau qui caractérise la classe de sol type 2 afin de déterminer le paramètre S. 
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Apres avoir calculé le produit aS , on effectue une interpolation linéaire pour obtenir les 
coefficients de réduction de G en tenant les écarts types, On obtient donc les résultats 


Rapport d’accélération 
as après l’interpolation |réducteur 


suivants : 

















Tableau 18: Les rapports d'accélération et le coefficient réducteur de calcul 


PD Vitesse de b nede Coefficient | Contrainte 
; densité du sol| . . cisaillement ; 4 
Les sites cisaillement Vs | Gsmax G; de poisson | admissible 
p (kg/m3) retenu G 
(m/s) (Mpa) y Oso (MPa) 


ICI 





, ,3 


Tableau 19 Les caractéristiques mécaniques des différentes catégories de site 
Remarque : 


Les contraintes admissibles du sol (déduites de DTU 13.12 Règles pour le calcul des 
fondations superficielles) seront utilisées dans le pré-dimensionnement des fondations 
pour chaque catégorie du sol. 


IV. Les études de cas à étudier : 
Les études de cas, qu’on va étudier, sont des systèmes portiques en béton armé qu’on va 
changer à chaque fois le nombre d’étages et les travées. 


1. 1% cas: 


C’est un portique simple d’un système plan de deux étages composé d’une seule travée : 
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K Ku 
Sans ISS Avec ISS 


Figure 35: Modélisation portique 2 étages avec et sans ISS 








C) 


Kg : Š Kh 





2. 2M cas: 


C’est un portique simple de quatre étages 








Ky 


Sans 155 Avec 155 


Ky 


Figure 36: Modélisation portique 4 étages avec et sans ISS 
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3. 3Me cas : 


C’est un portique simple de huit étages 





| K, Tk 


Figure 37: Modélisation portique 8 étages avec et sans ISS 
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4, 4M cas: 


C’est un portique à 3 travées de huit étages 









Sans 155 Avec 155 


Figure 38: Modélisation portique 8 étages et 4 travées avec et sans ISS 


V. Méthode de modélisation sur le logiciel ETABS : 


1. Présentation du logiciel ETABS : 


ETABS est un module de SAP2000, C’est un logiciel de génie civil qui permet la 
modélisation, l'analyse structurale et la conception des bâtiments. ETABS offre des outils 
de modélisation et de visualisation à base d'objets 3D inégalés, une puissance analytique 
linéaire et non linéaire claire, des capacités de conception sophistiquées et complètes 
pour une large gamme de matériaux et des écrans graphiques, des rapports, Et des dessins 
schématiques qui permettent aux utilisateurs de déchiffrer rapidement et facilement et 
facilement leurs résultats d'analyse et de conception. 
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2. Méthode de modélisation de PISS sur le logiciel: 


La modélisation de la structure avec ISS s’effectue en utilisant l’option « Link/Support 


properties » 


Define | Draw Select Assign Analyze 


[E Material Properties... 


B sms O A 


LEE Spring Properties b 


Diaphragms... 

Pier Labels... 
5pandrel Labels... 
Group Definitions... 
Section Cuts... 
Functions 


Generalized Displacements.. 


Display 
| 3-d Plá elf 





Detailing 
gep 35M É 


Frame Sections.. 


Design Options 


Tendon Sections... 


Slab Sections.. 
Deck Sections.. 


Wall Sections... 
Reinforcing Bar Sizes... 
Link/Support Properties... 


Frame/Wall Nonlinear Hinges... 


Panel Zone... 


Figure 39: étape 1 pour la modélisation sur le logiciel ETABS 


Puis, on clique sur l’option * Add a new property” afin d’insérer les valeurs de rigidités et 
les amortissements correspondants a ce “Link” après avoir bien sûr fixé la direction U3 et 
les rotations R1 et R3 vu que nos cas sont des portiques plans 


kii Linear Link/Support Directional Properties 





Link/Support Property Name 
ferme cent @) Stiffness ls Uncoupled 
u1 u2 
Directional Control 442463,996  538606,522 
Direction Fixed 


W] va E] Yes 

7] u2 E] Les rigidités 
u3 
R1 


R2 
v O) Damping ls Uncoupled 


] R3 | u1 u2 


1114,733 1235,899 
Shear Distance 


u2 0 


with respect to J-End of the 
link object 





Stiffness Values Used For All Load Cases 


Damping Values Used For All Load Cases 


Note: Distance is measured @S a m O rtisse m e nts 


Stiffness ls Coupled 


R1 R2 
Fixed 159165 3828 


Damping ls Coupled 


u3 R1 R2 
Fixed Fixed 663,4029 











Figure 40: étape 2 pour la modélisation sur le logiciel ETABS 


Ensuite, on définit un appui élastique ponctuel dans l’option « Spring properties » 
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Define | Draw Select Assign Analyze 
EN Material Properties... 


section Properties 


Spring Properties | 11 Point Springs. 


Diaphragms.. $3 Line Springs... 


Area Springs. 
Fier Labels... = rea prin gs 


Spandrel Labels... E Soil Profiles... 


E L, Isolated Column Footings... 
Group Definitions... 





Figure 41: étape 3 pour la modélisation sur le logiciel ETABS 


Puis, on associe le « link » défini précédemment à cet appui élastique en cliquant sur 
« Add » dans l’onglet du « Single Joint Links at point » 


pda Point Spring Property Data 


General Data 
Property Name 
Display Color 


Property Notes _ Modify/Show Notes... 


Spring Stifíness Options 


ía) User Specified Link Properties 


Simple Spring Stifness in Global Directions 
Translation X kN/m 
Translation Y kN/m 
Translation Z kN/m 
Rotation about X-Ads kN-m/rad 
Rotation about Y-fxis kN-m/rad 
Rotation about Z-Axis kN-m/rad 


Single Joint Links at Point 


Link Property dal Direction Axis 2 Angle 
| ferme cent g | +7 0 


N 


C'est le Link défini précédemment 





Figure 42: étape 4 pour la modélisation sur le logiciel ETABS 
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Remarque : 


On peut définir directement les rigidités des ressorts dans l’option « spring properties » 
sans les associer à un « Link » mais dans ce cas lá les amortissements ne seront pas 
tenus en compte. 


Tout d’abord, avant d’entamer les calculs des raideurs, on doit pré-dimensionner les 
semelles pour chaque cas grâce aux descentes de charges obtenues sur le logiciel 
« ETABS ». 


Les valeurs de la raideur globale de ressort dans les déférentes directions, pour chaque 
méthode, et pour chaque catégorie de site sont explicitées dans les tableaux suivants : 


VI. Résultats des raideurs et des amortissements 


Pour chaque cas, nous avons calculé les raideurs et les amortissements en adoptant les 
trois méthodes (NEWMARK-ROSENBLUETH , DELEUZE et le guide du SETRA) 


Le tableau 20 représente les dimensions des semelles retenues pour le calcul des 


Effort de Dimension de la semelle retenue a=b (m) 


compression P 
(KN) sol rocheux | sol meuble 


| 


Tableau 20: les dimensions des semelles et les efforts de compressions 


raideurs : 













ler cas : 
Méthode de NEWMARK-ROSENBLUETH : 


© Sites | K:(kN/m) | K,(kN/m) | Ko(kN.m/rad) | Ci(kNs/m) | C,(kNs/m) | Co (kN s/rad) 


1400000 1642666,667 182933,3333 923,186982 1161,343607 365,123951 
Sol rocheux 


84500 102142,8571 39232,14286 297,0560547 418,826534 95,22670187 
Sol meuble 
10080 13200 9720 125,459643 1946774683 36,53416586 
Sol mou 





Tableau 21: Résultats des raideurs et amortissements selon NEWMARK-ROSENBLUETH 


pour le 1er cas 
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Méthode de DELEUZE : 


Kn (kN/m) K, (kN/m) | Ko (kN.m/rad) | Cn (kN s /m ) | C, (kN s /m ) | Co (kN s /rad) 


1337237,279 1564653,662 157261,052 697,8703634 | 742,1474368 | 233,8446197 
Sol rocheux 


80233,25776 99842,83457 33576,51179 180,5704464 | 211,3108735 108,2291663 
Sol meuble 
9445,741301 12045,9407 8207,627604 | 24,61650748 | 1992552681 | 54,55720321 


Tableau 22: Résultats des raideurs et amortissements selon DELEUZE pour le 1er cas 





Méthode de guide SETRA : 


Ka (KN/M) K, (Nm) | Ko (Nam/rad) 


1516541,601 1685046,223 181379,5125 
Sol rocheux 
89289,13409 104778,0655 38898,90824 
Sol meuble 
10262,31158 13540,55001 9637,439113 
Sol mou 


Tableau 23: Résultats des raideurs et amortissements selon le guide du SETRA pour le 
1er cas 








Note : ; 
Les formules d’amortissements données dans le guide du SETRA ne semblent pas 
valides, les résultats étant aberrants. Utilisations des formules d’amortissement à 
éviter. Et par conséquent, nous allons comparer que les amortissements des 
méthodes de NEWMARK-ROSENBLUETH et DELEUZE. 


On compare la variation des raideurs et d’amortissement à l’aide des courbes qui 
traduisent cette variation en fonction de la méthode adoptée et le type du sol choisi : 
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RAIDEUR HORIZONTALE 


—— methode Newmark-Rosenblueth 8 methode de deleuze 


Methode du guide SETRA 


Leo D 


L 


— -g 
SOL MEUBLE SOL MOU 


RAIDEUR DE BALANCEMENT 


— methode Newmark-Rosenblueth 8 methode de deleuze 


Methode du guide SETRA 


SOL MEUBLE SOL MOU 


AMORTISSMENT HORIZONTAL 


—— methode Newmark-Rosentblueth =g methode de deleuze 


SOL ROCHEUX SOL MEUBLE SOL MOU 


AMORTISSEMENT DE BALANCEMENT 


—+= methode Newmark-Rosenblueth =g Methode de deleuze 


SOL MEUBLE SOL MOU 


2016-2017. 


RAIDEUR VERTICALE 


— methode Newmark-Rosenblueth 8 méthode de deleuze 


Methode du guide SETRA 


—— 
SOL MEUBLE SOL MOU 


—— méthode Newmark-Rosenblueth 


=J méthode de Deleuze 


—#— Méthode du guide SETRA 


AMORTISSEMENT VERTICAL 


—— methode Newmark-Rosentblueth =g methodededeeuze 


SOL MEUBLE SOL MOLI 


== méthode Newmark-Rosenblueth 


=f methode de Deleuze 





Figure 44: Comparaison les amortissements pour le 1er cas 
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Remarque : 


L”allure de la courbe des raideurs est presque la méme pour toutes autres cas, alors que 
pour les courbes des amortissements, 1l y a un changement de l’allure surtout pour les 
amortissements de balancement. 


by] 2ème cas : 
Méthode de NEWMARK-ROSENBLUETH : 


| Sites | K,(kN/m) K, (kN/m) Ko (KN.m/rad) | Cr (kN s/m ) | C, (kN s/m) | Co (kN s /rad) 
2000000 2346666,667 533333,3333 1639,657554 | 2106,04648 617,305719 
Sol rocheux 


117000 141428,5714 104142,8571 526,2542205 | 756,000764 158,075815 
Sol meuble 
14560 19066,66667 29293,33333 239,0291828 | 380,941588 62,7282045 


Tableau 24: Résultats des raideurs et amortissements selon NEWMARK-ROSENBLUETH 
pour le 2ème cas 





Méthode de DELEUZE : 


| Sites K, (kN/m) K, (kN/m) Ko (kN.m/rad) Ch (kNs/m) | CG (kNs/m) | Co (kN s /rad) 
1910096,786 2233480,795 458547,8667 1173,215149 | 1252,15734 564,701183 
Sol rocheux 
110955,7579 137177,6329 89248,07365 294.2686425 | 338,985633 249 428588 
Sol meuble 
13603,00967 18388,96073 24886,81969 48.05903102 | 259837573 133,108879 


Tableau 25: Résultats des raideurs et amortissements selon DELEUZE pour le 2ème cas 





Méthode de guide SETRA : 


Sites K, (kN/m) K, (kN/m) Ko (kN.m/rad) 
2166488,001 2407208,89 528803,2435 
Sol rocheux 


123631,1087 145077,3215 103258,2762 
Sol meuble 
14823,33895 19558,57223 29044,51815 
Sol mou 


Tableau 26: Résultats des raideurs selon le guide du SETRA pour le 2ème cas 
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AMORTISSMENT HORIZONTAL AMORTISSEMENT VERTICAL 


— méthode Newmark-Rosenblueth méthode de deleuze —— méthode New mark-Rosenblueth —S- méthode de deleuze 


SOL MEUBLE 


AMORTISSEMENT DE BALANCEMENT 


—+— méthode Newmark-Rosentblueth —— methode de deleuze 


K —— méthode Newmark-Rosenblueth 
: x méthode de Deleuze 


SOL ROCHEUX SOL MEUBLE SOL MOU 





Figure 46: Comparaison des amortissements pour le 2ème cas 


W] 3ème cas: 
Méthode de NEWMARK-ROSENBLUETH : 


K,(kN/m) | K,(kN/m) | Ko(kN.m/rad) | C,(kNs/m) | C,(kKNs/m) | Co (kNs/rad) 


2800000 3285333,333 1463466,667 2875,932188 | 3763,29421 1033,27804 

Sol rocheux 
169000 204285,7143 313857,1429 988,1253543 | 1458,66011 270,730476 

Sol meuble 
20720 27133,33333 84421,66667 450,4056917 | 734,588705 107,245196 


Tableau 27: Résultats des raideurs et amortissements selon NEWMARK-ROSENBLUETH 
pour le 3eme cas 





Méthode de DELEUZE : 


K, (kN/m) K, (kN/m) Ko (kN.m/rad) [C.(KNs/m)|C,(kKNs/m)| Co (kNs/rad) 


2673800,556 3124515,152 1258417,303 1952,777106 | 2091,53328 1322,96526 
Sol rocheux 


160109,5585 196928,9033 269251,125 491,8196078 | 559462698 612,471 777 
Sol meuble 
e 19389,61803 26748,74482 72025,07469 110,1094266 | 300,235696 318,956212 


Tableau 28: Résultats des raideurs et amortissements selon DELEUZE pour le 3ème cas 
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Méthode de guide SETRA : 


Kn (kN/m) K, (kN/m) Ko (kN.m/rad) 


3033083,201 3370092,446 1451036,1 
Sol rocheux 


178578,2682 209556,131 311191,2659 
Sol meuble 
21094,75159 27833,35279 83704,59592 
Sol mou 


Tableau 29: Résultats des raideurs selon le guide du SETRA pour le 3ème cas 





W] 4éme cas: 


RAIDEUR HORIZONTALE RAIDEUR VERTICALE 


—— méthode Newmark-Rosenblueth 4 méthode de deleuze —— méthode Newmark-Rosenblueth —m— méthode de deleuze 


Methode du guide SETRA Methode du guide SETRA 


mi "j 
m 1 i 


SOL ROCHEUX SOL MEUBLE SOL MOU SOL ROCHEUX SOL MEUBLE 


RAIDEUR DE BALANCEMENT 


—+— methode New mark-Rosenblueth 8 methode de deleuze 


—— méthode Newmark-Rosenblueth 


Methode du guide SETRA 


=J méthode de Deleuze 


Meéthode du guide SETRA 


SOL MEUBLE SOL MOL 





Figure 48: Comparaison des raideurs selon les 3 méthodes pour le 3ème cas 
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AMORTISSMENT HORIZONTAL AMORTISSMENT HORIZONTAL 


—— methode New mark-Rosenblueth —s— methode de deleuze —— méthode New mark-Rosenblueth —s— methode de deleuze 


SOL ROCHEUX SOL MEUBLE SOL MOU SOL ROCHEUX SOL MEUBLE 


AMORTISSEMENT DE BALANCEMENT 


—+— methode New mark-Rosenblueth =g methodededeeuze 


== méthode Newmark-Rosenblueth 


=f- methode de Deleuze 


SOL ROCHEUX SOL MEUBLE 





Figure 49: Comparaison des amortissements pour le 3ème cas 


W] 4éme cas: 
Méthode de NEWMARK-ROSENBLUETH : 


K,(kN/m) | K,(kN/m) | Ko(kN.m/rad) | C:(kNs/m) | C,(kNs/m) | Co (kNs/rad) 


4000000 4693333,333 4266666,667 5167,998568 6924,54606 1755,55371 
Sol rocheux 
240500 290714,2857 904517,8571 1835,640824 2776,71625 471,473538 
Sol meuble 
29120 38133,33333 234346,6667 838,8050613 1396,18985 192,098122 
Sol mou 


Tableau 30: Résultats des raideurs et amortissements pour les fondations centrales 
selon NEWMARK-ROSENBLUETH pour le 4ème cas 
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K,(kN/m) | K,(kN/m) | Ko(kN.m/rad) Ci (kN s /m ) C, (kN s/m) | Co (kN s /rad) 


2800000 3285333,333 1463466,667 2849,790556 3735,10097 | 1019,85999 
ol rocheux 


Sol meuble 169000 204285,7143 313857,1429 981,7292254 1451,77973 | 267482493 
Sol meuble 
| Sotmou 20160 26400 77760 430,8551745 700,277141 | 103,608006 
Sol mou 


Tableau 31: Résultats des raideurs et amortissements pour les fondations latérales 
selon NEWMARK-ROSENBLUETH pour le 4ème cas 





Méthode de DELEUZE : 


K,(kN/m) | K,(kN/m) | Ko(kN.m/rad) | Ci(kNs/m) | C,(kNs/m) |Co(kNs/rad) 


3820003,74 | 4465617,469 3668571,255 3301,009578 3528,07008 3184,5234 

Sol rocheux 
228033,0624 | 281648,5028 775305,3198 838,5070052 963,636251 1465,63455 

Sol meuble 
27182,26836 | 36293,00966 198866,3335 122,1774247 794,465619 751,014183 


Tableau 32: Résultats des raideurs et amortissements pour les fondations centrales 
selon DELEUZE pour le 4ème cas 





K,(kN/m) | K,(kN/m) | Ko(kN.m/rad) | Ci(kNs/m) | C,(kNs/m) |Co(kNs/rad) 


2673701,226 | 3123824,438 1258464,792 1923,877464 2062,72152 1305,44191 

Sol rocheux 
160082,7595 | 196727,7101 269297,9761 483,9371083 549,341962 604.366951 

Sol meuble 
18861,27356 | 25969,16277 66257,85834 98,74016089 332,639663 302,229098 


Tableau 33: Résultats des raideurs et amortissements pour les fondations latérales 
selon NEWMARK-ROSENBLUETH pour le 4ème cas 
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Méthode de guide SETRA : 


K, (KN/m) K, (kN/m) Ko (kN.m/rad) 


4332976,002 4814417,78 4230425,948 
Sol rocheux 
254130,6124 298214,4942 896834,9563 
Sol meuble 
29646,67791 39117,14446 232356,1452 
Sol mou 


Tableau 34: Résultats des raideurs pour les fondations centrales selon le guide du 
SETRA pour le 4ème cas 





KEN Ko Nan rad 


3033083,201 3370092,446 1451036, 1 
Sol rocheux 





178578,2682 209556,131 311191,2659 
Sol meuble 
olmo 20524,62317 27081,10001 77099,5129 
Sol mou 


Tableau 35: Résultats des raideurs les fondations latérales selon le guide du SETRA 
pour le 4eme cas 
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RAIDEUR DE BALANCEMENT 


—+— methode Newmark-Rosenblueth 8 méthode de deleuze 


Méthode du guide SETRA 


SOL MEUBLE 


SOL MOU 


2016-2017. 


RAIDEUR VERTICALE 


—— méthode New mark-Rosenblueth 8 méthode de deleuze 


Methode du guide SETRA 


ms 


i 


SOL MEUBLE 


—@— méthode Newmark-Rosenblueth 
=f méthode de Deleuze 


—à— Méthode du guide SETRA 


RAIDEUR VERTICALE 


= Méthode Newmark-Rosenblueth —- méthode de deleuze 


Methode du guide SETRA 


SOL MEUBLE 


—@— méthode Newmark-Rosenblueth 
= méthode de Deleuze 


fr Méthode du guide SETRA 





Figure 51: Comparaison des raideurs latérales selon les 3 méthodes pour le 4eme cas 
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Figure 52: Comparaison des amortissements centraux pour le 4ème cas 
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Figure 53: Comparaison des amortissements latéraux pour le 4ème cas 
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Discussion : 


% 


ee 


% 


ee 


On remarque que les valeurs des raideurs sont trés proches entre eux. 

Les raideurs obtenues de la méthode de DELEUZE sont légérement petits par 
rapport aux autres méthodes. 

On remarque également une diminution des amortissements en fonction de la 
déformabilité du sol et une augmentation de ceux-c1 avec l’élancement du 
bâtiment. 

Les amortissements obtenus de la méthode de DELEUZE sont inférieurs de ceux 
de NEWMARK-ROSENBLUETH, à part l’amortissement de balancement où on 
constate un dépassement de la méthode de Deleuze qui commence à partir du 1% 
cas pour un sol mou. 


Conclusion : 


On peut dire que la méthode de DELEUZE prend en compte plusieurs 
paramètres, notamment la fréquence qui n’est pas prise en considération par les 
autres méthodes. Cependant, cette méthode exige beaucoup de calcul et 
d’interpolation puisqu'elle dépend des transmittances, et par conséquent, notre 
étude sera restreinte sur la méthode de NEWMARK-ROSENBLUETH car c’est la 
méthode la plus utilisée vu sa simplicité, et ne nécessite pas de calculer le rayon 
équivalent dans le cas des semelles rectangulaires. 

La présence de l’amortissement a pour effet d’augmenter la compression du sol. 
En effet, l’amortissement permet d’entraver tout mouvement du déplacement et 
de basculement du sol, et donc le sol doit être comprimé pour limiter son 
déplacement, c’est pour cela que l’amortissement diminue en fonction de la 
déformabilité du sol. 

Plus la hauteur du bâtiment est grand, plus on a un effort de compression au 
niveau du sol important, ce qui implique l’allongement de l’amortissement avec 
l’élancement du bâtiment. 

l’effet du non linéarité des sols dans l’Interaction Sol-Structure augmentait 
essentiellement dissipation (amortissement) de l’énergie. 
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Chapitre VIII : Résultats et interprétation sur les cas étudiées 


Ensuite, on va entamer le calcul de la période fondamentale, le déplacement en téte, 
déplacement inter-étages et l’effort tranchant à la base sur le logiciel à l’aide des valeurs 
des raideurs et d’amortissements calculés précédemment, on obtient pour chaque cas les 
résultats suivants : 


I. Résultats obtenus de la méthode NEWMARK- 
ROSENBLUETH : 


Y] 1% cas: 













Sol rocheux | Sol meuble | _sormou _[sotrocheux | Sot meuble] sormon 


Tableau 36: Les résultats trouvés pour le 1èr cas 





La variation de période fondamentale pour les deux cas (avec et sans ISS) est représentée 
par la courbe suivante : 


PÉRIODES AVEC ET SANS ISS 
(NEWMARK-ROSENBLUETH) 


——aveciSS sans ISS 
0,7 
0,6 — 
0,5 
0,4 
0,3 


0,2 
0,1 


SOL ROCHEUX SOL MEUBLE SOL MOU 


Figure 54: Évaluation de la 1ère période en fonction de site avec et sans ISS pour le 1er 
cas 
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Concernant les déplacements latéraux pour chaque type du sol : 


la 


Sol Rocheux Sol meuble Sol mou 


=p sans |SS =f avec ISS 
Figure 55: variation des déplacements en fonction de site avec et sans ISS pour le 1er cas 


L’effort tranchant à la base est représenté par le graphe suivant : 


Effort tranchant a la base (Newmark-Rosenblueth) 


120 
100 


80 


60 
40 
20 
0 Dh | E — 


sol rocheux sol meuble sol mou 





Esans ISS E avec ISS 


Figure 56: Effort tranchant à la base en fonction du site avec et sans ISS pour le 1èr cas 
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Type du site es Sans ISS | Avec ISS epiac ment Vérification 
étages admissible (mm) 
étage 1 | 2619 | 2,636 15 | OK | 
roeas [étagel [ 2619 | 2636 | 
DC 


O 
OK 
m s o 





Tableau 37: Vérification des déplacements inter-étages avec et sans ISS pour le 1er cas 


Puis, on va effectuer une comparaison en étudiant la variation de la période 


fondamentale, le déplacement en tête et l’effort tranchant à la base pour chaque site : 


Comparaison des périodes 
1,004395604| 1,08571429| 1,4175824 
Comparaison des déplacements 


1,005838469 | 1,19412175| 1,4293203 
Comparaison des efforts tranchants 
1,002125783| 1,09172857| 1,4068732 


Tableau 38: Comparaison des déplacements et les efforts tranchants pour chaque 
catégorie du site pour le 1er cas 





Commentaires : 


+ 
ey 


Dans le cas sans ISS, la période propre reste fixe peu importe le site, et dans le cas 
avec ISS il y a une augmentation de la période en fonction du site, on constate pour 
cette cas une augmentation de 0,4% pour le site 1, 8,5% pour le site 2 et 41,76% pour 
le site 3. 

Pour le déplacement, il y a un allongement pour les deux cas : avec et sans ISS. De 
plus, cette augmentation est notable dans le cas avec ISS surtout pour un site trés 
meuble (mou). La comparaison entre le site sans ISS et avec ISS montre que le 
déplacement en téte augmente de 0,5% pour le site 1, 19,4% pour le site 2 et 42,93% 
pour le site 3. 

Le déplacement inter-étage qui le déplacement relatif entre deux étages successifs est 


inférieur au déplacement admissible ( daam = 0,01 x = avec h=3m :la hauteur de 


l’étage et K=2 : c’est le coefficient de comportement ) pour les deux cas et pour touts 
les sites. 
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Y 


* On constate que plus le sol est déformable, plus l’effort tranchant à la base augmente 
dans les deux cas : sans ISS et avec ISS. Cependant, il y a une réduction de cet effort 
des sites en ISS par rapport aux sites avec ISS de 0,2% pour le site 1, 9,2% pour le 
site 2, et 40,7% pour le site 3. 


W] 2ème cas: 



















Tableau 39: Les résultats trouvés pour le 2ème cas 








PÉRIODES AVEC ET SANS ISS 
(NEWMARK-ROSENBLUETH) 


—@— avec ISS sans ISS 


1,4 
1,2 

1 
0,8 
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0,4 
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Figure 57 : Évaluation de la 1ère période en fonction de site avec et sans ISS pour le 
2ème cas 
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Figure 58: Variation des déplacements en fonction de site avec et sans ISS pour le 
2ème cas 
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Figure 59: Effort tranchant à la base en fonction du site avec et sans ISS pour le 2ème 
cas 
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Tableau 40: Vérification des déplacements inter-étages avec et sans ISS pour le 2ème 
cas 


Comparaison des périodes 
1,004716981| 1,16037736| 1,49646226 
Comparaison des déplacements 


1,007613877| 1,20273209| 1,65093369 
Comparaison des efforts tranchants 
1,003226055 |  1,10989574| 1,40975016 


Tableau 41: Comparaison des déplacements et les efforts tranchants pour chaque 
catégorie du site pour le 2ème cas 





Commentaires : 


Y 


% Dans le cas sans ISS, la période propre reste fixe peu importe le site, et dans le cas 
avec ISS il y a une augmentation de la période en fonction du site, on constate pour ce 
cas une augmentation de 0,4% pour le site 1, 16,03% pour le site 2, et 49,64% pour le 
site 3. 

* Pour le déplacement, il y a un allongement pour les deux cas : avec et sans ISS. De 

plus, cette augmentation est notable dans le cas avec ISS surtout pour un site mou La 

comparaison entre le site sans ISS et avec ISS montre que le déplacement en téte 

augmente de 0,7% pour le site 1, 20% pour le site 2, et 65,1% pour le site 3. 
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Y 


* Le déplacement inter-étage qui le déplacement relatif entre deux étages successifs est 
pros 2 Ed h 
inférieur au déplacement admissible ( daam = 0,01 X z avec h=3m : la hauteur de 
l’étage et K=2 : c’est le coefficient de comportement) pour les deux cas et pour touts 
les sites. 

% On constate que plus le sol est déformable, plus l’effort tranchant à la base augmente 

dans les deux cas : sans ISS et avec ISS. Cependant, il y a une réduction de cet effort 


des sites en ISS par rapport aux sites avec ISS de 0,3% pour le site 1, 11% pour le site 
2, et 41% pour le site 3. 


W] 3ème cas: 













Sans ISS Avec ISS 


Sol mou Sol mou 
Effort tranchant à la base F (kN) 160,303 113,3556 











Tableau 42: Les résultats trouvés pour le 3ème cas 


PÉRIODES AVEC ET SANS ISS 
(NEWMARK-ROSENBLUETH) 


—@— aveciSS =f sans ISS 
2,5 
1,5 


0,5 


SOL ROCHEUX SOL MEUBLE SOL MOU 


Figure 60: Évaluation de la 1ère période en fonction de site avec et sans ISS pour le 
3ème cas 
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Figure 61: Variation des déplacements en fonction de site avec et sans ISS pour le 3ème 
cas 


Effort tranchant à la base (Newmark-Rosenblueth) 


180 
160 
140 
120 
100 
80 
60 
40 
20 


0 Y 1 : 
sol rocheux sol meuble sol mou 





MW sans ISS E avec ISS 


Figure 62: Effort tranchant à la base en fonction du site avec et sans ISS pour le 3ème 
cas 
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ica an Les EEI Sans ISS | ISS | Avec ISS Lopes Vérification 
ica MET ON SEE 


EA 
DS E MES E AI ES 


Z e E 
a EE 088 a 161 EI 
E E A e 
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Tableau 43: LL des déplacements inter-étages avec et sans ISS pour le 3*"* 
cas 
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Comparaison des périodes 
1,006535948 | 1,17290553|  1,64052288 
Comparaison des déplacements 


1,009992608 1,33029828| 1,9074596 
Comparaison des efforts tranchants 
1,00436437 1,17583911 1,41416039 


Tableau 44: Comparaison des déplacements et les efforts tranchants pour chaque 
catégorie du site pour le 3ème cas 





Commentaires : 


% Dans le cas sans ISS, la période propre reste fixe peu importe le site, et dans le cas 
avec ISS il y a une augmentation de la période en fonction du site, on constate pour ce 
cas une augmentation de 0,6% pour le site 1, 17,3% pour le site 2, et 64,05% pour le 
site 3. 

* Pour le déplacement, il y a un allongement pour les deux cas : avec et sans ISS. De 
plus, cette augmentation est notable dans le cas avec ISS surtout pour un site mou. La 
comparaison entre le site sans ISS et avec ISS montre que le déplacement en tête 
augmente de 1% pour le site 1, 33,02% pour le site 2, et 90,74% pour le site 3. 

* Le déplacement inter-étage est admissible pour le cas sans ISS, mais n’est pas 
admissible pour le cas avec ISS pour le site 4. 

% On constate que plus le sol est déformable, plus l’effort tranchant à la base augmente 

dans les deux cas : sans ISS et avec ISS. Cependant, il y a une réduction de cet effort 

des sites en ISS par rapport aux sites avec ISS de 0,4% pour le site 1, 17,6% pour le 

site 2, et 41,41% pour le site 3. 


W] 4ème cas : 










Sans ISS Avec ISS 
Sol meuble Sol meuble 


Période fondamentale T(s) 1,928 1,928 1,928 1,932 2,087 2,319 
Déplacement absolu (mm) 45,999 64,383 82,835 46,104 72,088 105,2 
Effort tranchant à la base F (KN) 241,3164 337,843 434,3696 240,93 74 319,2654 | 371,8707 


Tableau 45: Les résultats trouvés pour le 4ème cas 
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PÉRIODES AVEC ET SANS ISS 
(NEWMARK-ROSENBLUETH) 


—@— aveciSS sans ISS 


2,5 


2 m— y . 


1,5 


0,5 
SOL ROCHEUX SOL MEUBLE SOL MOU 


Figure 63: Évaluation de la 1ère période en fonction de site avec et sans ISS pour le 


s 
4eme cas 
o 86385883 0585383232238 _= 8 8 2 3 Y 
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Figure 64: Variation des déplacements en fonction de site avec et sans ISS pour le 
4ème cas 
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Effort tranchant a la base (Newmark-Rosenblueth) 
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Figure 65: Effort tranchant à la base en de“ du site avec et sans ISS pour le 4ème 


an 
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| étage 5 | 8,415 | 11,503 | 15 | OK | 





Tableau 46: Vérification des déplacements inter-étages avec et sans ISS pour le 4ème 


cas 


Comparaison des périodes 
1,002074689 | 1,08246888| 1,20280083 
Comparaison des déplacements 


1,002282658 | 1,11967445| 1,26999457 
Comparaison des efforts tranchants 
1,001573023| 1,05818858| 1,16806621 


Tableau 47: Comparaison des déplacements et les efforts tranchants pour chaque 
catégorie du site pour le 4eme cas 





Commentaires : 


Dans le cas sans ISS, la période propre reste fixe peu importe le site, et dans le cas 
avec ISS il y a une augmentation de la période en fonction du site, on constate pour 
cette cas une augmentation de 0,2% pour le site 1, 8,2% pour le site 2, et 20,3% pour 
le site 3. 

Pour le déplacement, 1l y a un allongement pour les deux cas : avec et sans ISS. De 
plus, cette augmentation est notable dans le cas avec ISS surtout pour un site mou. La 
comparaison entre le site sans ISS et avec ISS montre que le déplacement en téte 
augmente de 0,23% pour le site 1, 12% pour le site 2, et 27% pour le site 3. 

Le déplacement inter-étage est admissible pour le cas sans ISS, mais 1l n’est pas 
admissible pour le cas avec ISS pour le site 4.De plus, le nombre d’étages dont le 
déplacement inter-étages n’est pas admissible a diminué. 

On constate que plus le sol est déformable, plus l’effort tranchant à la base augmente 
dans les deux cas : sans ISS et avec ISS. Cependant, il y a une réduction de cet effort 
des sites en ISS par rapport aux sites avec ISS de 0,16% pour le site 1, 5,81% pour le 
site 2, et 16,8% pour le site 3. 
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Conclusion : 


Il y a un allongement de la période fondamentale de la structure lorsque l'effet de 
PISS est pris en compte par rapport a celle obtenue lorsque la structure est considérée 
encastrée à sa base. Ce qui implique théoriquement une atténuation des efforts 
sismiques et une variation plus ou moins de l’accélération de la structure. 

Dans le cas avec ISS, le déplacement de la structure augmente en fonction des 
sites sous l’effet de l’interaction sol structure. On observe une augmentation de celle- 
ci plus le sol est déformable. En effet, une augmentation de T correspond aussi à une 
augmentation des déplacements de la structure par rapport au sol. On note également 
que la flexibilité et la souplesse de la structure permettent d’augmenter ce 
déplacement. 

Effet de l’inertie de la structure : les efforts tranchants à la base des structures fixées à 
leur base sont plus importants que ceux des structures dont la base est élastique. Dans 
ce cas précis, on constate que l’effort tranchant à la base augmente quand le nombre 
d’étage augmente. En effet, une partie de l’énergie est absorbée par le sol, et 
l’amortissement du système couple bâtiment-sol réduit considérablement les 
accélérations dans la superstructure. L'hypothèse de l’encastrement parfait de la 
structure dans le sol n’est vérifiée que lorsque le sol est comporte une grande rigidité 
et la réponse dans ce cas ne que du chargement extérieur statique ou dynamique 

Les déplacements inter-étages ne sont pas admissibles pour le 3°" et le 4° cas pour 
un sol mou dans le cas avec ISS, alors que dans le cas sans ISS ces déplacements sont 
admissibles, ce qui conduit à un sous-dimensionnement si on considere un bâtiment 
élancé qui est encastré dans un site mou, et par conséquent, des dommages 
imprévisibles vont être générés dans la structure. 

Le fait d'augmenter le nombre des travées (4% cas) permet de diminuer les effets de 
PISS sur notre structure puisqu'il y a une diminution de la variation des paramètres 
par rapport au 3%" cas. Donc, on peut dire que le nombre de travées du portique 
permet de rigidifier notre structure. 

Le fait de prendre en considération l’ISS n’a pas un effet notable sur les sites rocheux 
et ferme, vu la petite variation des paramètres qu’on a calculé auparavant, et donc, on 
peut négliger les effets de l’ISS pour ces deux types de site et considérer notre 
structure comme parfaitement encastrée dans le sol. Alors que dans les sites meubles 
et mous on doit tenir en compte les effets de l’ISS surtout pour les structures élancées 
ou des structures qui sont caractérisées par leurs masses, dans ce cas lá, on doit 
prendre en considération la flexibilité du sol, sinon on aura un mauvais 


dimensionnement des paramétres de notre structure. 
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IL —P”effet de amortissement sur la structure : 
On va étudier l’effet de l’amortissement sur l’effort tranchant à la base et les 
déplacements pour les 4 cas, pour un sol mou du modèle avec ISS, vu que c’est le cas le 
plus défavorable, on obtient les résultats suivants : 


Avec Sans Avec Sans 
amortissement | amortissement | amortissement | amortissement 


RU 814 nu ,002 44 [RO Es 441 












Déplacement en tête (mm) 
Effort tranchant à la base (kN) D or En 1,1828 Maca 101 Misal 104 


Avec Sans Avec Sans 
amortissement | amortissement | amortissement | amortissement 


RL 68 144, Mia E ” EN 82 











Déplacement en tête (mm) 
Effort tranchant à la base (kN) | un 0 Mill a 1,8707 Md 


Tableau 48: L'effet de l'amortissement sur le déplacement en tête et l'effort tranchant 
à la base pour les 4 cas 


Discussion : 


L’amortissement radial est modélisé dans le matériau par un amortisseur visqueux 
(énergie dissipée proportionnelle à la vitesse du mouvement). La présence de cet 
amortissement visqueux a pour effet d'augmenter la résistance en compression du sol, 
entravant tout mouvement de basculement à mesure que sa valeur augmente. À la vue de 
ces résultats, nous pouvons conclure que l’amortissement n’a pas une influence notoire 
sur le comportement globale de la structure, dans cette comparaison, on constate que plus 
l’amortissement est élevé plus les efforts tranchants à la base et les déplacements 
diminuent, mais d’une manière qui n’est pas significative. En effet, l’influence 
l’amortissement devient prépondérante avec la plasticité du sol, ce quí n’est pas notre cas 
(sol élastique linéaire). 


IT. Prise en compte de PISS en utilisant la méthode forfaitaire: 


Lorsqu'un bâtiment comporte une infrastructure, on peut évaluer de façon approximative 
l’effet de l’interaction sol-structure en définissant une « hauteur de dimensionnement » H 
supérieure à Ho hauteur de la structure au-dessus du sol. H définit en fait le niveau 
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supposé de l’encastrement parfait ; 1l inclut une partie du sous-sol et traduit la flexibilité 
de ce sous-sol. La période fondamentale T est plus élevée, les déplacements sont plus 
grands que si on considère un encastrement au niveau du sol. 

Selon le règlement PS92 pour un bâtiment qui comporte une infrastructure 1l faut prendre 
en considération la hauteur de celle-ci dans les calculs pour avoir l'effet de PISS, cette 
hauteur varie selon le type du sol: 


e Pourun sol de très bonne résistance (sol rocheux) : 
H — Ho 


e Pour un sol de résistance mécanique moyen (sol meuble) : 


Hı 
Ge 


e Pour un sol de résistance mécanique faible (sol mou) : 
H — Ho + Hı 
Avec 


Ho: Hauteur de la structure au-dessus du sol. 


H, :Hauteur du sous-sol. 

















Figure 66: Hauteur de dimensionnement d'après PS92 


N.B: 

Cette méthode permet seulement d’avoir la période propre en tenant compte de 
Pinteraction sol-structure, elle ne permet pas d’évaluer l’effort tranchant à la base car 
avec l'augmentation de hauteur la masse augmente est par suite l’effort tranchant a la 
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base augmente ce qui est contradictoire avec ce qu’on a trouvé dans la modélisation par 
des ressorts. 


Or si on veut le calculer on peut adopter la méthode statique équivalente à partir des 
périodes données par la méthode forfaitaire. 


1. Périodes : 


On a effectué les calculs pour le 3°" et le 4ème cas en changeant le nombre des sous- 
sols et leurs hauteurs, on obtient les périodes suivantes pour chaque type de sol: 


3ème cas (R+8 avec un sous-sol, une 
seule travée) 


Méthode ressort Méthode forfaitaire 


SETORES) [Soussolde3mf 197 | 2761 | L906 [2135 | 


Tableau 49: Comparaison des périodes selon la méthode de ressort et la méthode 
forfaitaire en fonction des sous-sols ( 3** cas) 



















Commentaire 1: 

Lorsqu’on est en présence d’un sol meuble, la période obtenue de la méthode forfaitaire 
s’approche radicalement de la modélisation par des ressorts surtout dans le cas de deux 
sous-sols. Alors que dans le sol mou, il y a une différence notable des périodes : la 
méthode forfaitaire tend a s’approcher à la méthode de ressorts avec l augmentation des 
sous-sols et leurs hauteurs. 





4®™e cas (R+8 avec un sous-sol, trois 
travées) 


Méthode ressort Méthode forfaitaire 


périodes(s) | Un soussotdeám | 2087 [259 | 205 | 241. 


Tableau 50: Comparaison des périodes selon la méthode de ressort et la méthode 
forfaitaire en fonction des sous-sols (4*"* cas) 














Commentaire 2: 

Pour cette cas, on constate que les périodes sont très proches pour les deux sites et les 
deux méthodes, surtout pour le cas d’un seul sous-sol de hauteur de 4 m où on peut 
approximer la modélisation par des ressorts avec la méthode forfaitaire. 
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2. Déplacements : 


3e cas 8 étages 


Tableau 51: Comparaison des déplacements selon la méthode des ressorts et la 
méthode forfaitaire en fonction des sous-sols (3°"° cas) 













Commentaire 1: 


Comme pour les périodes, les déplacements en tête de la structure avec la méthode 
forfaitaire s’approchent mieux à ceux de la méthode des ressorts avec l’augmentation 
de hauteurs et des nombres d’étages, on observe qu’ils s’approchent radicalement pour 
une structure de deux sous-sols de 4m. 






4éme cas 8 étages/ 3 travées 


[meuble | mou | meuble | mou | meuble | mou 


Tableau 52: Comparaison des déplacements selon la méthode des ressorts et la 
méthode forfaitaire en fonction des sous-sols (4°"° cas) 







Commentaire 2: 


De même que pour les périodes de ce cas, on constate que les déplacements sont très 
proches pour les deux sites et les deux méthodes, surtout pour le cas d’un seul sous-sol de 
hauteur de 4 m. 


Discussion: 


D’après ce qui précède, on peut dire que la méthode forfaitaire donne des résultats 
comparables avec la modélisation par des ressorts, cette méthode est valable pour les 
bâtiments réguliers qui ont des sous-sols. En effet, les résultats obtenus dans la 4°" cas 
sont très proches entre les deux méthodes en raison que c’est une structure qui s’approche 
de la réalité, au contraire, du 3% cas qui s’éloigne de la réalité, ce qui interprète la 
différence obtenue entre les deux méthodes surtout pour un sol de mauvaise résistance 
(sol mou). 
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Conclusion : 


Les conditions géologiques ainsi que les propriétés des sols doivent faire l’objet d’une 
attention spéciale car elles peuvent être sources de dégâts irrémédiables. On remarque 
que la modélisation d’une structure avec tenir compte l’effet ISS sont les plus fiables car 
leur dépendance des caractéristiques des sols est trop importante. 
Alors que la modélisation d’une structure encastrée (sans tenir compte l’effet ISS) n’est 
pas à recommander car elle donne des résultats trop conservatifs. 
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Chapitre IX: Fondations profondes 


Concernant les fondations profondes il n’existe pas beaucoup de formules permettant de 
calculer la raideur d’un pieu d’une manière directe et simple, pour notre étude on va 
s’intéresser à la méthode du Fascicule 62 Titre V et l'EUROCODE 8 pour évaluer la 
raideur en tête du pieu. 


I. Détermination des raideurs selon Fascicule 62 Titre V : 
Lorsque l'on a affaire à un sol hétérogène composé de plusieurs couches avec des 
caractéristiques variées la modélisation la plus adaptée pour le comportement d'un pieu 
est le modèle de Winkler (Figure 62 ci-dessous). Ce modèle consiste à modéliser le pieu 
comme une poutre verticale reposant sur une série de ressorts indépendants horizontaux 
(latéral). Le Fascicule 62 Titre V Fondations profondes donne des formules permettant de 


déterminer les raideurs des ressorts à partir des résultats des essais au pressiomètre. 


Mouvement du sol en champ 


Ñ Mouvement du pieu 
libre P 


Ressort latéral 


Couchej 


MY 





Figure 67: Modèle dynamique de Winkler pour la modélisation de 
l'interaction sol-pieu 


1. Raideur horizontale :( L’annexe C5 du Fascicule 62 Titre V) 


La méthode décrite dans L’annexe CS, article 3 du Fascicule 62 Titre V Fondations 
profondes : loi d’interaction vis-à-vis des sollicitations de courte durée d’application 
considère que le sol exerce en chaque section de l'élément une réaction perpendiculaire à 
l'axe de celui-ci, qui est fonction du déplacement transversal relatif de la section. Cette 
réaction se compose de : 
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e Pressions frontales, modélisées par une pression uniforme s'exerçant sur la 
plus grande largeur de l'élément perpendiculairement au sens du déplacement 


e Efforts de frottement transversaux sexercant sur les parties du périmètre 
paralléles au sens du déplacement 


Dans le cas de pieux circulaires on considére que la réaction se compose uniquement de 
pressions frontales. 


La loi de mobilisation de la réaction frontale est donnée en fonction du déplacement ò du 
pieu est définie par : 


e Un segment de droite passant par l’origine et de pente Kf 


e Unpalierr: 





Figure 68 : Réaction frontale en fonction 
du déplacement ó 


L’évaluation de Kret de rı sont obtenus à l’aide des formules suivantes : 


e Pour Kf qui est le module linéique de mobilisation de la pression frontale pour un 
élément de fondation profonde: 


12E 
kem _ Lorsque B > Bo (sollicitations de courte durée D'application) 
Ptas j 
3 ' Bo 
12Em O , a 
Kf —;—— Lorsque B < Bo (sollicitations de courte durée D'application) 


=X(2.65)%+a 


e Pour le palier rl : 
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Avec: 


Em le module préssiométrique de Ménard. 

B est la largeur de l’élément perpendiculaire au sens du déplacement. 
Bo est une largeur de référence prise égal à 0,60m. 

a est le coefficient rhéologique. 

Pf* est la pression de fluage nette. 


KKK SK S$ 


Sous-consolidé altéré 
el remanié ou lâche 





Tableau 53: les valeurs de a. et le rapport Em/P, selon le fascicule 62 


Grâce à ces formules 1l est possible de déterminer les raideurs de chaque couche de sol. 


NB : Cette formule est valable tant que l’on n 'atteint pas le palier rl sinon, il y aura des 


déformations plastiques du sol. 


2. Raideur verticale : 


L'annexe G4 du Fascicule 62 Titre V donne la méthode suivante : 


Pour entamer les calculs 1l faut admettre que l’élément de la fondation se comporte de 
façon élastique et linéaire, et donc la rigidité axiale de l’élément est prise égale : 


R; = 2c (Sollicitations de courte durée D application) 


ci 


Avec 


Qc: Charge de fluage du pieu. 


S ci : représente l’enfoncement provoqué par la charge Qg supposé de courte durée d’application. 


2.1.Calcul de Penfoncement Spi : 
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L’enfoncement provoqué par la charge Qc est définie par la relation suivante: 


QcPL 


€; = 
li E;S 





B 
Sci = Tog + À 


: représente le raccourcissement instantané sous la charge Qe de la partie de 


Pélément non comprise dans la hauteur d’encastrement et supposée non soumise à 


frottements. 


Ir 
> 
1 





Figure 69: représentation de l'enfoncement des pieux 
sous une charge Qc 


D, : Longueur libre. 


Avec : 


D:Hauteur du pieu. 


D, :Hauteur d’encastrement équivalent. 





1 D 
D, = — x | P,"(z)dz 
Le d 
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E; : Module d’élasticité du béton instantané (= 32 000 MPa pour un B25) 


E, = 110007 |f., 


S: Section du pieu. 


2.2.Calcul des limites de fluage: 


Les charges de fluage d'un élément de fondation profonde pour les éléments «mis en 
ceuvre sans refoulement du sol sont évaluées a partir de Qpu et de Q.,, par les relations 
suivantes : 


- En compression : Qe = 0.5Qpu + 0.7Q 
- Entraction Qie = 0.704 


Qpu : Charge limite de rupture en pointe ; 


Qu : Charge limite de rupture au frottement ; 
2.2.1. Calcul de la contrainte de rupture sous pointe : 
La contrainte limite de pointe est donnée par : 
Apu = Kp X Pie 
Avec : 


Pie : pression limite nette équivalente. 
- kp: facteur de portance, fonction de la nature de la formation concernée et du 
mode de mise en œuvre de l'élément de fondation. 
La charge limite de pointe sera donc : Qp= A.Qpu . Avec A: la section de pointe du 
pieu. 


La pression limite nette équivalente p,, est calculée par l’expression : 


1 D+3xa 
pi == *) pt (2) x dz 
° b+3xa Jos, 


Avec : 


- a: max (B/2 et 0,50 m). pour notre cas c’est 0.5 
- b:min (a, h). égale a “a” et donc de valeur 0.5m 
- h: ancrage dans la couche où se situe la pointe du pieu qui vaut : 
o Pour une seule couche de sol h = 6B avec un minimum de 3m, 
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o Pour un sol multicouche pour lequel la contrainte effective cv” due au 
poids des terrains au-dessus de la couche d’ancrage est au moins égale à 
100kPa (environ 7 à 10m de terrain) h = 3B avec un minimum de 1.5m ; 
c’est la règle classique des 3 diamètres d’ancrage du pieu dans la couche 


résistante 





Figure 70: illustration des paramètres de calcul 


p*i(z): est obtenu en joignant linéairement par des segments de droites les 
différentes valeurs mesurées. 


Le choix du facteur de portance se fait en se basant sur les deux tableaux suivants : 


Valeurs du facteur de portance k, 


Éléments Éléments 


Nature des terrains mis en œuvre sans mis en œuvre avec 
refoulement du sol refoulement du sol 


Craies 


C 
Mames. Marno-calosires 18 26 | 
Roches altérées (1) IAB 1,8 à 3,2 


(1) La valeur de k, pour ces formations est prise égale à celle de la forma- 
tion meuble du tableau à laquelle le matériau concerné s'apparente le 
plus. Le calcul proposé ici n'affranchit, en aucun cas, d'avoir recours aux 
méthodes spécifiques de la mécanique des roches. 





Tableau 54: Valeurs du facteur de portance Kp 
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Definition des catégories conventionnelles de sols 


ver | 
Argiles etlimons mous <07 | <30 
Ardiles etlimons formes 12820 | 3060 
Ailes très fermes à dures | >25 | >60 
[e es E 
|Moyennement compacts 10820 | BOB ABO 
Compas E | »200 


Argiles. Limons 


Sables. Graves 


OA Mol eo E 
Craies OB aes AT 

~ B Compas o o 
(1) L'appellation de roches altérées ou fragmentées peut regrouper des matériaux calcaires, schisteux ou d'origine granitique. S'il est difficile parfois de fixer des 


IEC 
Compactos IE IA 
| Tendres tsaso 
rames MOTS a 
Atérées 25880 OOO 

Roches (1) | 
Fragmentées RES EE 
limites précises avec les sols meubles qui constituent leur phase finale d'évolution, on réservera toutefois cette classification aux matériaux qui présentent des 

modules pressiométriques supérieurs à 50 à 80 MPa. 





Tableau 55: Définition des catégories conventionnelles de sols 


2.2.2. Calcul du frottement latéral unitaire limite : 


L’effort total mobilisable par frottement latéral sur la hauteur h concernée du fût du pieu 
se calcule par : 


h 
Qu = PX | as(z) x dz 
0 


P : est le périmètre du pieu. 


Le frottement latéral unitaire limite est donné par l’abaque suivant : 








2 


E RS SSD | 
PRELELEELE EESTE: TANE 
3 4 5 


př (MPa) 





Figure 71: L'abaque de qs en fonction de pl* 
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Le choix de la courbe se fait sur le tableau suivant : 


Choix des courbes pour le calcul du frottement latéral unitaire q, 
(principaux types de pieux utilisés en ouvrage d'art) 


Nature du sol 


Type de pieu Argiles. Limons Sables. Graves 


aam à Jajasujela an] a 


Foré del 


CI 
Foré à la boue BE ser fa 07,112) | Gp, Gp | Or | Os | Q4.Qsí1) | Os | Osas | Os 


(tube récupéré) 


Jena perdu) 


Réalésage et rainurage en fin de forage. 

Pieux de grande longueur (supérieure à 30 m). 

Forage à sec, tube non louvoyé. 

Dans le cas des craies, le frottement latéral peut être très faible pour certains types de pieux. Il convient d'effectuer une étude spécifique dans chaque cas. 
Sans tubage ni virole foncés perdus (parois rugueuses). 

Injection sélective et répétitive à faible débit. 

Injection sélective et répétitive à faible débit et traitement préalable des massifs fissurés ou fracturés avec obturation des cavités. 





Tableau 56 : Chois des courbes pour le calcul du frottement latéral unitaire 


II.  Eurocodes : 


La rigidité est définie comme la force ou le moment qu'il faut appliquer à la tête du pieu 
pour produire un déplacement ou une rotation unitaire suivant la même direction et nulle 
suivant les autres directions. L’annexe C de l’Eurocode 8-5 donne les rigidités statiques 
équivalentes en tête du pieu pour trois types de modèle de sol : 


e un sol dont le module d'Young varie linéairement avec la profondeur 


DL 
= Esx= 


e un sol dont le module d'Young varie avec la racine carrée de la profondeur 


E=E;x | 
= Es F 


e un sol dont le module d'Young reste constant avec la profondeur 


E = Es 
Les formules de rigidité horizontale Kyy, rigidité à la flexion Kym et rigidité de 
couplage Kym = Kmx sont données dans le tableau suivant : 
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Modèle de sol K y (MPa.m) K um(MPa.m?) © Kyy = Kmu(MPa.m”) 
Z E E E 
E = Eç x — d x Es x 0.6 x (035 d? x Es x 0.14 x (008 | =d? x Es x 0.17 x ($6 
d Es Es Es 


E E E 
E = Eş X d x Es X 0.79 x e d? x Es x 0.15 x TI —d? x Es X 0.24 x Ti 


AIN 


E E E 
E=Es d x Es x 1.08 x (021 d? x E; x 0.16 x (20758 | =d? x Eç x 0.22 x (250 
Es Es Es 


Tableau 57: Les formules de rigidité horizontale Kpn, rigidité à la flexion Kmm et rigidité 
de couplage Kum=Kmn 





Avec: 


e Ele module d'Young du modèle de sol, égal à 3G 

e Ez le module d'Young du matériau constitutif du pieu 

e Es le module d'Young du sol à une profondeur égale au diamètre du pieu 
e dle diamètre du pieu. 

e z la profondeur. 


La matrice de raideur en tête de pieu s'écrit avec les notations précédentes : 


ME Kan Km >i 
M Kun Kmum Ile 


X et ® désignent le déplacement horizontal et la rotation en tête de pieu 
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HI. Exemple de calcul des rigidités d’un pieu : 


Les données de calcul : Pieu de diamètre 1.2m avec un module d'Young de 32000MPa et 
d’une longueur de 37,5 


| Profondeur Em 


Remblai 
a ë 
75 

[9 | 0204 | 1.55 [7598039216] 05 | 

a 

Sable grossier légèrement vaseux propre —— A oai 
A 


Tableau 58: Les données de calcul du pieu 
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1. Raideurs selon le Fascicule 62 Titre Y : 


Tout d’abord on va calculer les Kẹ 


Profondeur (m) Ke 


74,72356776 
65,84376173 
6 | 13,14211293 

-7,5 13,36410808 
[9 | 9,158298368 
9,867327919 
8,567440408 
31,5011119 
8,803783592 
12,25413232 


94,83422135 
145,4692296 
151,7323239 
350,9696278 
136,7836176 
166,3265155 
191,6152362 
378,3658866 
383,053215 

373,1268012 
450,5291939 
445,565987 

541,0486332 
477,70866 

389,2572235 


Tableau 59: calcul de rigidités différées en fonction de la profondeur 





Et pour le calcul de la rigidité équivalente, on va modéliser le pieu sur ETABS comme 
une poutre verticale, puis on va affecter sur chaque 1.5m de profondeur un ressort de 
rigidité Kf, et par la suite on va appliquer une force unitaire (un moment) à la tête du pieu 
ce qui va engendrer un déplacement (une rotation). 
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3000 


AER EEE A | 





Figure 72: Modélisation du pieu sur le logiciel ETABS 
Les rigidités équivalentes sont données par les formules suivantes : 


1.1.Rigidité horizontale : 


A _F 
A 


Avec: F = 3000KN Et A= 23,218mm 


AN: Ky = 129,21MN /m 


1.2.Rigidité de rotation : 


a M 
RD 
Avec: M = 10000KN.m Et € = 0.012495rad 
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AN: Kg = 800,32MN.m/rad 


1.3.Calcul de la rigidité verticale : 


R; = ze (Sollicitations de courte durée D application) 
ci 


Avec: Qc = 0.5 * Qu + 0.7 * Qs 


On obtient les résultats illustrés dans le tableau suivant : 


Profondeur (m) 


IS 168 


A 


Tableau 60: Les pressions limites mesurées par l'essai pressiométrique 
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Alors 


1 39.3 
Pie = 77 X f pi (z) x dz = 3542,5KPa 
2.4 36.9 


D’après le tableau 50 la classe du sol : sable, grave type C compacts. 


Y Facteur de portance k, : 
D’après le tableau 49 : 


k, = 1.2 (Pieu foré donc pas de refoulement du sol) 


Donc 


qpu = 1.2 X 3542.5 = 4251KPa 
Qpu = A X q, = T X 0.6" x 4251 = 4807,77KN 


v Calcul Pression limite nette équivalente p;, sur la hauteur du pieu: 
1 D 
Pie = D X f pi (z) x dz = 1799,49KPa 
0 


D’après tableau 50 la classe du sol : sable, grave type B moyennement compact. 


Donc d’après le tableau 51 pour un pieu foré à la boue dans un sable, gravier de type B 
on est dans la courbe Q1 (pieu de grande longueur supérieure à 30m). 


Et par la suite d’après la courbe de la figure 65 : 
qs = 0.04MPa 
Alors 
Qu = T X 1.2 x 37.5 X 0.04 x 1000 = 5654,87KN 
Ainsi: 


Qc = 0.5 x 4807,77 + 0.7 x 5654,87 = 6362.29KN 


B _B , xD, 
100 E¡xS 
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Avec D, = D — D, 





Donc 
D; = 37.5 — 19,05 = 18.45m 
E; = 32000MPa 


AN: 


1.2 6362.29 x 18.45 


-E p CMA AAA 01524342 
100 * 32000 x 1000 x x x 0.62 de 


Sci 


Alors la rigidité verticale du pieu 


K, = Ca — 417,37 MN 
VES.  0,02724364  ” ¡00 


Le tableau suivant regroupe les résultats obtenus: 


Rigidité horizontale 129,21 MN/m 


Rigidité de rotation 800,32MN.m/rad 
Rigidité verticale 417,37 MN/m 


Tableau 61: tableau récapitulatif des raideurs selon le fascicule 62 





2. Raideurs selon L'EUROCODE 8: 
Données du calcul 


Pieu de diamètre 1.2m avec un module d'Young de 32000MPa et d’une longueur de 
37.5m. 


Donc 


E, = 32000MPa Et d = 1.2m 
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Or Eoea = m, Module cedométrique 


Avec E y : Module pressiométrique 
2v? ; 
E = Ejea(i — 1): Module d'Young 


Avec v : Coefficient de poisson du sol (on prend v = 0.33 pour cet exemple) 


Profondeur (m) | EM (MPa) | «| EouMPa) | EMPa) 
| 29% | 0.67 | 4417910448 | 296 
9 f 15 |] 05 | a | 2077 


Tableau 62: Calcul des modules d'Young a partir des modules pressiométriques 





Et par suite : E = E(1,2m) = 15.16 x = = 12.48MPa (On considère une variation linéaire 
entre 0 et 1.5m) 


D’autre part on observe que ce cas s’approche plus au 2eme modèle dont le module d'Young 
varie avec la racine carrée de la profondeur : 


E= Esx | 
= Es a 
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Résultats trouvés : 


Modele de sol Ky(MPa.m) Kum(MPa.m3) | Kyy = Kuym(MPa.m?) 
E=Esx E 106,54 1363.73 -276,39 


Tableau 63: tableau récapitulatif des raideurs selon l'EC8 


Tableau récapitulatif: 


Fascicule 62 Titre V 


Rigidité horizontale | Rigidité de rotation | Rigidité verticale 
(MN/m) (MN .m/rad) (MN/m) 


129,21 800,32 417,37 


Rigidité horizontale | Rigidité à la flexion Rigidité de 
(MN/m) (MPa.m') couplage (MPa.m?) 
106,54 1363,73 276,39 


Tableau 64: comparaison de raideurs trouvées selon le fascicule 62 et l'EC8 





Conclusion: 


Contrairement aux fondations superficielles, il ne semble pas exister de formules générales 
qui permettent de déterminer simplement et rapidement la raideur d’un pieu. 


Concernant les formules de 1*EC8, on a pu constater d’une part qu'il était très rare d’avoir un 
sol dont les caractéristiques sont constantes en fonction de la profondeur. Il est également rare 
que celles-ci augmentent linéairement en fonction de la profondeur. Il est ainsi très difficile de 
classer les sols dans l’une des trois catégories donnés par l’Eurocode or les différentes 
formules donnent des raideurs qui peuvent varier jusqu’à 30%. D'autre part, L’EC&8 ne prend 
pas en compte la longueur du pieu. 


Ainsi, il semble plus simple et plus réaliste de modéliser le pieu dans un logiciel avec des 
ressorts. En effet, ces derniers permettent de prendre en compte chaque couche de sol avec ses 
caractéristiques propres (les Kf calculés par Fascicule 62 Titre V) ce qui n’est pas le cas dans 
1?EC8.Une fois le pieu modélisé il suffit de lui appliquer un effort et à partir des déplacements 
qui en résultent de calculer les raideurs en tête des pieux pour revenir à la même approche 
traitée auparavant. (Chaque pieu est modélisé dans la structure par des ressorts suivant chaque 
degré de liberté et dont leurs raideurs sont ceux calculées en tête du pieu) 
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Chapitre X : Application sur une structure 
tridimensionnelle 


I. Présentation de la structure : 


e La structure qu’on va étudier est un bâtiment d’habitation (classe II) de huit 
niveaux (R+7) et située dans la zone de Rabat (zone II). 

e C’est un bâtiment régulier qui est composée d’une ossature en portique en béton 
armé qui repose sur un sol de résistance faible (sol mou). 

e La dalle considérée est une dalle pleine d’épaisseur 25 cm. 


e Les dimensions des poteaux après un pré dimensionnement: 


La section des poteaux (cm x cm) 


Base 65x65 


Tableau 65: Prédimensionnement des différents poteaux 





e La dimension de la section des poutres : 40x80 (cm x cm) 
e La portée des poutres est de 8 m, la hauteur des poteaux est de 3m. 
e Les charges surfaciques appliquées sont : 


v Charge permanente G=250 kg/m? 
Y Charge d’exploitation Q=250 kg/m? 


6 E37 bos Bo Bo 
> PE bas os E; À 5 = 
RA À ` : 
À 7) - À < 
> ES $ SS < EA 
E Mo Bo RL M: À À À | / j A 
RSR SLI 
AS >, 2 Es 
: E A a LA 
y | | PPI l 
B E R; E; de | j 
y SY / 
<< Sy 
j 
| | ; 
Sac $, +, t, $, | 


E X M R; Ebo Lo 


Figure 73: Modélisation de la structure sur ETABS 
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IL. Calcul des raideurs de ressort : 


Les résultats obtenus sont illustrés dans le tableau suivant : 























Raideur Raideur de Amortissement 
































P Raideur vertical Raideur de Amortissement Amortissement Amortissement 
horizontale balancement ; 7 ; de balancement . 
Ka Ko torsion Ko horizontal Ch vertical Cy Co de torsion Co 
(KN/m) (KN/rad) (KN.s/m) (KN.s/m) (KN.s/rad) (KN.s/rad) 


(KN/rad) 





Tableau 66: Calcul des raideurs et d'amortissements 


En raison de la symétrie de la structure,et aprés avoir prédimensionné les semelles de la 
structure on aura trois catégories de ressort (notés : mou 1 , mou 2 et mou 3) qu’on va 
calculer leurs raideurs à l’aide de la méthode de NEWMARK-ROSENBLUETH. 





T mou! mou2 mou2 mou2 at 
F mou2 mou3 mou3 mou3 mou2 
T mou2 mou3 mou3 mou3 mou2 
pi 


+) mou2 mou3 mou3 mou3 mou2 


F mou2 mou3 mou3 mou3 mou2 
Y 


— mou? X mou2 mou2 mou2 mou 


Figure 74: Disposition des différents ressorts 
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IHI. Comparaison des résultats obtenus sur le logiciel : 


1. Période et effort tranchant : 


Les périodes et l’effort tranchant à la base : 









| Périod e(s) | Di CU rt tranc ha: at( KN) | 
Tableau 67: Calcul des périodes propres et l'effort tranchant pour les 2 cas selon les 
deux directions x et y 


Il y a une augmentation de la période 30% suivant x et de 26 % suivant y. Cependant, il y 
a une réduction de l’effort tranchant à la base de 22% suivant x et de 19% suivant y. 


2. Les déplacements : 


2.1. Les déplacements latéraux : 


| Déplacements latéraux À (mm) 
Suivant y 
Sans ISS 
61,736 l i 
56,94 : : 
48,305 | | 
36,744 : : 


22,282 ; 
16,668 | 
10,253 : 


Tableau 68: Les déplacements latéraux pour les 2 cas 


30,34 43,794 | 30,029 | 43,204 





suivant x suivant y 
== sans [SS =f avec ISS 
Figure 75: Représentation des 


déplacements latéraux pour les 2 cas (avec 
et sans ISS) 


selon les directions x et y 
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Il y a un allongement du déplacement, notamment le déplacement en tête qui atteint 28% 
suivant x et 25% suivant y. 


e Déplacement admissible A est égal : 
Aadm= 0.004 X 26 x 1000 = 104mm 


Tous les déplacements sont admissibles pour les deux cas. Il est à noter qu'il y a un 
allongement notable pour le cas avec ISS. 


2.2. Les déplacements inter-étages : 








Déplacement inter-étages d (mm) 
Entre 5 ot 6% étage 
Entre 39% et 45% étage 
Entre 1% et2 étage 
Tableau 69: Comparaison des déplacements inter-étages pour les 2 cas selon les 
directions x et y 









Déplacement admissible inter-étage d est : 
daam = 0.005 x 3 x 1000 = 15mm 


Tous les déplacements inter-étages sont admissibles. De même, 1l y a un allongement des 
déplacements inter-étages en raison de l’augmentation des déplacements latéraux. 


3. Dimensionnement des éléments porteurs de la structure : 


Dans cette partie on va dimensionner les poutres et les poteaux pour les deux cas ( avec et 
sans ISS) afin de comparer les sections obtenus et déduire l’effet de PISS sur les 
sections d’acier. 


3.1.Les numérotation adoptées : 


Les numérotations adoptées par le logiciel pour les éléments porteurs : 
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vV Les poteaux : 








i ru k os 4. Hu 

5 Hilo >; La mo M: 

Te Mo Mo M M: 

3) Fa AE M; AT M7 

Fo E; E; mt M: 
Y 

a s di Ms E7 Es E, 


Figure 76 : Numérotation des poteaux 


Y” Les poutres : 








B26 





B29 


Figure 77: Numérotation des poutres 
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> Les résultats de dimensionnement obtenus sont rassemblés dans l’annexe 1 et 2. 
3.2.Le poids d’acier obtenu : 


On va s'intéresser seulement sur la variation de la section d’acier tout en fixant la 
section du béton. 


3.2.1. Les poteaux : 


Le poids d’acier obtenu pour les poteaux pour les deux cas avec et sans ISS : 


Poids d’acier Poids d’acier 
(cm2) (cm2) (m) (sans ISS) (kg) (avec ISS )(kg) 
582,4 


Tableau 70: Le poids d'acier pour les poteaux obtenu pour les deux cas 





















3.2.2. Les poutres : 


Le poids d’acier obtenu pour les poutres pour les deux cas avec et sans ISS : 


As (avec Poids d’acier Poids d’acier 
(cm2) ISS) (cm2) (m) (sans ISS) (kg) (avec ISS)(kg) 
Etage7 | 2258 | 2674 | 8 | 1418024 | 1679272 
| Etageó | 218 | 2513 | 8 | 136904 | 1578164 
| Etages | 2076 | 2465 | 8 | 1303778 | 158648 7 
| Etage4 | 2093 | 2342 | 8 | 1314404 | 1470778 | 
| Etage3 | 2011 | 222,8 | 8 | 1262008 | 13991848 7 
btage2 | 1932 | 2087 | s | im200 | neos | 
| Etagel | 187 | 192 | 8 | 117436 | 1209768 7 
Base | 1922 | 2095 | 8 | 1207016 | 131566 


Tableau 71: Le poids d'acier pour les poutres pour les deux cas 
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Discussion : 


On constate qu'il y a une diminution de volume d’acier dans les poteaux et une 
augmentation dans les poutres si on prend en compte l’interaction sol-structure. 


On peut bien visualiser la répartition des sections d’acier d’une manière détaillée dans les 
annexes : il y a une augmentation de la section d’acier pour les poutres et les poteaux 
centraux, et une diminution de celle-ci pour les autres poteaux, on peut justifier cela par 
le fait que le système d’appui a changé et par conséquent la répartition globale et locale 
va être changée. 


Y] Quantité d’acier totale : 


e Poids volumique acier=7850kg/m? 
e Prix du kg d'acier=12Dhs 












| 3450437296 | 3422027224 __ 


Tableau 72: Comparaison des prix d'acier obtenus pour les deux cas 





> Iln’y a pas une différence notable pour les deux cas en termes de prix. 


L”Eurocode 8-5 impose l’étude des effets de l’interaction dynamique sol-structure dans 
les cas suivants: 


+ 


% Structures pour lesquelles les effets du P-d (2e ordre) jouent un rôle significatif, 
ce qui est le cas pour les structures élancées, pour lesquelles l’excentrement du 
poids des masses par rapport à l'axe verticale induit un moment fléchissant 
supplémentaire. 

% Structures avec fondations massives ou profondes, comme les piles de ponts, les 

silos ou tout ouvrage industriel équivalent. 

Structures hautes et élancées, comme les tours et les cheminées. 

e Structures supportées par des sols très mous, sol de classe Sl avec des vitesses 

des ondes de cisaillement de Vi max<100 m/s. 


% 


ee 


ve 


Pour les structures pour lesquelles la prise en compte de PISS n’est pas imposée par PEC 
8-5, la prise en compte de l’ISS permet en général une réduction des sollicitations par une 
dissipation au niveau du sol. 
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On peut classer notre structure dans la catégorie des ouvrages supportés par des sols très 
mous. Dans ce cas l’effet de l’ISS sera défavorable : 


On remarque globalement qu'il y a une augmentation des sections de l’acier, et par 
conséquent, la non prise en compte de PISS conduit a une sous-estimation des valeurs de 
la section d’acier et donc à un sous-dimensionnement de notre ouvrage. Ainsi, la 
répartition des efforts dans les membrures va être différente en raison du changement de 
système d’appui ce qui va affecter sur la sécurité du bâtiment. 
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Conclusion et recommandation 


L’effet de l'action du sol sur le comportement de la structure lors d'un séisme joue 
un róle primordial et ne peut étre négligé ni dissocié du calcul de la réaction totale de la 
construction face à un séisme. L’hypothèse de l’encastrement parfait de la structure dans 
le sol n’est vérifiée que lorsque le sol comporte une grande rigidité et la réponse dans ce 
cas ne que du chargement extérieur statique ou dynamique. Ainsi, L'objet de notre 
présent travail a consisté sur la compréhension des aspects les plus importants de l'ISS en 
considérant une étude générale sur des portiques plans implantés sur des différents sites. 
En effet, L’étude de l'influence de l'interaction sol structure sur les paramètres 
définissant la réponse sismique a montré que celle-ci peut allonger d’une manière 
significative les périodes et les déplacements avec la flexibilité du sol et de réduire les 
efforts tranchants à la base. 

Ensuite, pour la modélisation de l’effet de PISS pour les fondations profondes, nous 
avons utilisé la méthode définie par l’annexe C5 du fascicule 62 titre V et la méthode de 
l'EUROCODE 8-5 pour évaluer la rigidité du pieu. La méthode de l'EUROCODE 8-3 est 
approximative en raison qu’elle se base sur trois modèles et elle ne prend pas en 
considération les différentes couches au contraire de la méthode du fascicule 62 qui prend 
en compte chaque couche de sol avec leurs propres caractéristiques. 

Enfin, une structure tridimensionnelle a été considérée afin de bien visualiser l’effet de 
l'interaction sol structure sur un ouvrage réel reposant sur un sol de mauvaise résistance 
(sol mou). A partir des résultats trouvés, la non prise en compte de PISS conduit à un 
sous-dimensionnement de la structure ce qui va affecter la sécurité de l’ouvrage. Donc, à 
partir des conclusions ci-dessus et des résultats obtenus, 1l apparait que la non prise de 
PISS conduit à des erreurs importantes qui pourraient conduire à une ruine ou à une 
dégradation importante de la structure concernée. 


A partir des résultats obtenus, il apparait que la non prise de PISS conduit à des erreurs 
importantes qui pourraient conduire à une ruine ou à une dégradation importante de la 
structure concernée, de ce fait il est recommandé de tenir compte de PISS pour les 
structures reposant sur des sols de mauvaises résistances. Il est recommandé également 
que l'interaction sol structure ne doit pas être négligée car elle est liée 
au type de sol et de la structure. Il est donc important de mener une étude de sol détaillée 
afin de pouvoir donner un jugement correct sur les caractéristiques exactes du sol et donc 
sur le type de site de calcul à considérer. De plus, le fait de tenir compte L’ISS est plus 
sécuritaire et n’a pas un effet sur le budget du projet. 


Et finalement, ce travail n’est que le début vers d’autres études plus approfondies. Il 
serait très fructueux d'introduire d’autres complexités en tenant en compte par exemple la 
plasticité du sol qui affecte d’une manière sensible l’interaction sol-fondation-structure et 
qui va nous imposé à adopter un comportement non linéaire du sol, et par conséquent, la 
méthode des éléments finis est indispensable dans l’analyse dynamique de la structure. 
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Annexes 
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Annexe 1 : Dimensionnement des poteaux 


Sans ISS : 


NEN) |M(KN.m)|MKN.m)| Mmw | Asicm”) | As(em’) | Asítotale)(em”)_ 
SE A A A A A | me 
| C2-C12-C23-C25 |772,9259| 134,5086 | 162,9074 | 162,9074 | _6 | 6 | 
n Goo rss 1286260 | 1393240 | 1286260 | 25 | 25 | ss 
o | cece.corcas |402.7418| 203,942 | 1775609 | 203942 | 12 | 12 | A > 
oo| crc [390/2984 207,168 | 185,7421 2071168 | 123 | 123 | 246 
| C10-C13-C22-C5 |388,3032| 126,2675 | 234,6001 | 234,6001 | 146 | 146 | 292 | 
y Cimc24  |768,4428| 134,5173 | 117,8794 | 134,5173 | 37 | 37 | 74 | 
| C14-C17-C18-C21 |372,0314| 124,6813 | 245,4557 | 245,4557 | 156 | 156 | 312 | 
C15-C16-C19-C4 |761,0722| 128,2893 | 169,856 | 169,856 | 66 | 66 | 132 | 
| C1-C9-C26-C30 |446,5739| 223,3941 | 221,5857 | 223,3941 | 132 | 132 | 264 | 
C2-C12-C23-C25 |1515,042| 225,1341 | 234,6586 | 234,6586 | 146 | 146 | 292 | 
W | C3-C20 |1509987|217,1193 | 190,5827 | 217,1193 | 129 | 129 | 258 | 
q | C6-C8-C27-C29 [845,5087 | 271,4685 | 211,8511 | 271,4685 | 145 | 145 | 29 
oo| creas  |833,7674| 2742128 | 183,5735 | 2742128 | 107 | 147 | 294 
4> | C10-C13-C22-C5 |786,7941| 216,4862 | 283,5396 | 283,5396 | 157 | 157 | 314 
HI | C11-C24  |1527,419| 225,1405 | 189,1388 | 225,1405 | 137 | 137 | 274 | 
| C14-C17-C18-C21 |757,6447| 211,7732 | 293,9389 | 293,9389 | 16,7 | 167 | 334 | 
C15-C16-C19-C4 |1503,516| 216,8305 | 241,5603 | 241,5603 | 151 | 151 | 302 | 
| C1-C9-C26-C30 | 715,0422 | 272,4336 | 253,6499 | 272,4336 | 16,7 | 167 | 334 7 
C2-C12-C23-C25 |2284,179| 286,6608 | 281,3217 | 286,6608 | 248 | 248 | 496 | 
un | C3-C20 | 2254,78 | 278,8823 | 240,998 | 278,8823 | 23,9 | 239 | 478 | 
q | C6-C8-C27-C29 | 1314,01 | 314,1012 | 258,8576 | 314,1012 | 202 | 202 | 404 
oo| czes —|1295/013| 316092 | 233.1186 | 3160922 | 202 | 202 | aa 
4> | C10-C13-C22-C5 |1193,317| 278,1274 | 308,4099 | 308,4099 | 188 | 188 | 376 | 
HA | C11-C24  |2298,781|286,7134 | 239,6065 | 286,7134 | 249 | 249 | 498 | 
| C14-C17-C18-C21 |1134,099| 273,0121 | 317,4351 | 317,4351 | 19,8 | 198 | 396 | 
C15-C16-C19-C4 |2248,416| 278,6694 | 287,9304 | 287,9304 | 246 | 246 | 492 | 
C1-C9-C26-C30 |1012,336| 308,0067 | 331,2434 | 331,2434 | 95 | 95 | 19 | 
C2-C12-C23-C25 |3041,432| 364,4882 | 429,233 | 429,233 | 209 | 209 | 418 | 
st | C3-C20 |3001,764|346,1777 | 334,1953 | 346,1777 | 149 | 149 | 298 | 
D [cree [1804641] 3975071 | 3714571 | 3975071 | 12 | 12 | n 
oo| ceas [1784306 4009492 | 316,5762 | 400.9992 | 123 | 121 | 242 
C10-C13-C22-C5 |1606,887| 342,4662 | 452,2579 | 452,2579 | 145 | 145 | 29 | 
| C11-C24 |3064,486| 364,7898 | 330,7454 | 364,7898 | 16,6 | 166 | 332 | 
| C14-C17-C18-C21 |1525,256| 330,2276 | 471,9505 | 471,9505 | 154 | 154 | 308 | 
| C15-C16-C19-C4 |2987,914| 345,7832 | 445,5902 | 445,5902 | 216 | 216 | 432 | 


Etage 


Etage 
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C1-C9-C26-C30 
C2-C12-C23-C25 
C3-C20  |3744751| 381,3192 | 356,6278 | 381,3192 | 23 | 23 | 46 | 
C6-C8-C27-C29 404,4627 
C7-C28 406,3556 14,8 
C10-C13-C22-C5 431,8 15,1 
C11-C24 395,7891 24,6 24,6 49,2 
C14-C17-C18-C21 447,2962 15,6 15,6 31,2 


UT 
UT 


LO + + y N 
LO UT IS y 
D Ea 00 00 O0 
O 00 00 00 


C15-C16-C19-C4 381,081 451,5337 27,5 27,5 
C1-C9-C26-C30 382,346 4,4 
C2-C12-C23-C25 600,3579 
C3-C20 456,112 11,6 
C6-C8-C27-C29 495,2319 5,9 
C7-C28 457,809 4,4 
C10-C13-C22-C5 546,378 
C11-C24 65,379 12,8 


N 


`~ 
`~ 


N 
O 


E 


11,6 
5,9 
4,4 


23,2 
11,8 


`~ 


15,4 
25,6 
17,8 
41,6 


`~ 


D 


W 


12,8 
8,9 
20,8 


C14-C17-C18-C21 571,711 8,9 
C15-C16-C19-C4 626,7321 20,8 
C1-C9-C26-C30 412,018 4,4 
C2-C12-C23-C25 616,1186 26,5 
C3-C20 473,007 18,4 
C6-C8-C27-C29 533,169 9,5 
C7-C28 445,077 
C10-C13-C22-C5 560,352 9,3 
C11-C24 490,699 20,2 


à E 
WO D N D 
WO UT y OR 


26,5 
18,4 
9,5 


UT 
Uy 


36,8 


E 


`~ 


9,3 
20,2 
10,2 
27,1 


18,6 
40,4 
20,4 
54,2 


C14-C17-C18-C21 |2728,116| 435,319 585,8705 585,8705 10,2 
C15-C16-C19-C4 15199,959| 471,1839 | 642,6711 642,6711 27,1 


C1-C9-C26-C30 10,3 10,3 20,6 
C2-C12-C23-C25 21,5 21,5 3 
C3-C20 20,1 20,1 40,2 
C6-C8-C27-C29 | | 19,8 
C7-C28 9,9 9,9 19,8 
C10-C13-C22-C5 9,7 
C11-C24 852,141 | 901,9418 21,6 
C14-C17-C18-C21 9,6 
C15-C16-C19-C4 890,1324 20 | 2% | 40 _ 
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> 
< 
© 
C> 
pnd 
7 
Y 


C2-C12-C23-C25 | 752,8822 | 161,4297 | 162,9074 | 1629074 | 6 | 6 
C3-C20 | 730,9827 | 135,8459 | 119,3244 | 135,8459 | 4 | 4 7 
C6-C8-C27-C29 28,2 
c7.c28 | 400,031 | 240,2954 | 185,7421 | 240,2954 | 15 | 15 
C10-C13-C22-C5 


?) J 
C11-C24 752,5296 | 167,1197 | 117,8794 | 167,1197 
7 7 


2 
12 
0 
28,6 
13 
7 
18 
9 


m) 
C1-C9-C26-C30 218,312 169,012 171,9604 | 171,9604 2 
, 3 


C14-C17-C18-C21 15,4 15,4 30,8 
C15-C16-C19-C4 | 732,0838 | 134,0489 | 169,856 | 169,856 | 67 | 67 | 134 
C1-C9-C26-C30 25,4 
C2-C12-C23-C25 28,2 
C3-C20 23,6 


NEO [MARNE [ma ETES 
11 11 
15 15 
6, 6, 
i 6, 6, 


C6-C8-C27-C29 30,2 
C7-C28 31,6 
C10-C13-C22-C5 30,8 
C11-C24 29,6 
C14-C17-C18-C21 32,6 
C15-C16-C19-C4 29,2 
C1-C9-C26-C30 26,8 
C2-C12-C23-C25 4 
C3-C20 43,2 
C6-C8-C27-C29 38,8 
C7-C28 40,4 
C10-C13-C27-C5 38,6 
C11-C24 48,4 
C14-C17-C18-C21 39,6 
C15-C16-C19-C4 474 
C1-C9-C26-C30 | 1091,064 | 325,3756 | 331,2434 | 331,2434 | 9 | 9 
C2-C12-C23-C25 3 


C3-C20 32,4 
C6-C8-C27-C29 29,4 
C7-C28 31,6 
C10-C13-C27-C5 29,6 
C11-C24 41,2 
C14-C17-C18-C21 31,6 
C15-C16-C19-C4 41 
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N(KN) | My(KN.m) | M-(KN*m) M max Ası(cm?) | As2(cm?) | Asítotale)(cm?) 


C1-C9-C26-C30 15,4 
C2-C12-C23-C25 35,8 
C3-C20 45,6 

C6-C8-C27-C29 32,2 
C7-C28 34,2 

C10-C13-C22-C5 31,2 
C11-C24 52,4 

C14-C17-C18-C21 33,2 
C15-C16-C19-C4 51,4 
C1-C9-C26-C30 | 
C2-C12-C23-C25 34,4 
C3-C20 

7 7 

9 9 


1 1 


00 
00 


15,6 
8 


C6-C8-C27-C29 2792,928 | 534,8421 | 495,2319 | 534,8421 


C7-C28 2848,272 | 555,7162 | 408,8565 | 555,7162 | 9 | 9 | 
C10-C13-C22-C5 | 2633,195 | 493,1353 | 546,3786 | 546,3736 | 8 | 8 | 


A A 

4 4 
C11-C24 4246,077 | 617,9889 | 450,5667 | 617,9889 | 18,8 18,8 37,6 
C14-C17-C18-C21 | 2618,79 | 466,4217 | 571,7113 | 571,7113 18,4 
C15-C16-C19-C4 | 4155,214 | 518,2088 | 626,7321 | 626,7321 37,4 
C1-C9-C26-C30 | 2231,813 | 399,8397 | 412,0181 | 412,0181 
7 7 
7 7 


2 
1 
1 


00 
00 


7 7 
8 
8 8 
6 
2 2 
C2-C12-C23-C25 44,8 
C3-C20 39,8 
C6-C8-C27-C29 20,2 
C7-C28 23 
C10-C13-C27-C5 20,2 
C11-C24 48,8 
C14-C17-C18-C21 22,6 
C15-C16-C19-C4 47,8 
C1-C9-C26-C30 10,4 
C2-C12-C23-C25 
C3-C20 
C6-C8-C27-C29 10,4 
C7-C28 10,4 
C10-C13-C27-C5 10,4 
C11-C24 34 
? ; A ; 2 2 


C14-C17-C18-C21 | 3622,392 | 595,9702 | 578,8016 | 595,9702 10,4 
C15-C16-C19-C4 | 5489,102 | 868,8954 | 873,0791 | 873,0791 32,4 


32,4 
34 
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Annexe 2 : Dimensionnement des poutres 


Elément 
B1-B5-B21-B25 148,5979 | 4 
6,9 6,5 


pa 
ay B11-815 251,0413 241,3145 
QD B12-B14 182,4652 141,7289 
+ 162,6228 125,6568 
B26-B29-B46-B49 | 192,4568 167,0788 
B27-B28-B47-B48 | 149,0751 | 4 | 85,6236 
B30-B33-B42-B45 | 280,225 255,7128 
B31-B32-B43-B44 | 215,8898 146,3977 | 4 
B34-B37-B38-B41 | 286,0654 258,39 
B35-B36-B39-B40 | 221,2277 144,6305 
B1-B5-B21-B25 136,2147 132,9579 
B2-B4-B22-B24 118,5402 83,7648 
B3-B23 104,7835 73,4372 
B6-B10-B16-B20 | 220,7805 | 6 | 216,5441 
B7-B9-B17-B19 179,0616 144,1233 
© B8-B18 156,4002 126,1816 
x B11-B15 222,7009 219,7688 | 6 | 
00 B12-B14 184,0246 144,7071 
+ 159,6092 125,5402 
B26-B29-B46-B49 | 179,336 166,6476 


LN 
[8284822824 | 116,517 | 31 | 815309 | 3 
le 
L| s6-B10-B16-820 | 2223271 | 61 | 2185306 | 6 
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B7-B9-B17-B19 176,4647 4,8 142,1305 3,9 
__ 811815 | 2240057 | 61 | 2193568 | 6 

-B39- | | 3,8 
__82-84-822824 | 109,2506 | 3 | 790038 | 3 
__ 83-823 | 100,0877 | 3 | 744955 | 3 7 
__86810816820 | 180,2999 | 49 | 1849637 | 5 
m 
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B35-B36-B39-B40 | 210,2123 141,3109 | 38 


Etage 4 


Un | © | 


Etage 


151 


Ecole Mohammadia d'Ingénieurs- Projet de Fin d'Etude 2016-2017. 


B31-B32-B43-B44 193,7776 5,3 136,4101 3,7 
2 


ll Bess | 149,622 | 4 | 1309995 | 36 
ll Bees | 171815 | 4 | 1314251 | 36 
B13 148,7009 | a | 1293276 | 36 | 

D 

(ge) 

(aa) 


Etage 2 


Etage 1 


ll 83-823 | 967099 | 3 | 76152 | 3 7 
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Uy 


Ecole Mohammadia d'Ingénieurs- Projet de Fin d'Etude 2016-2017. 


B8-B18 150,8444 134,8308 
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